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PREFACE 


This is the final report on research carried out under the two 
year (Jan. 1983 to Dec. 1984) cooperative agreement NCC1-76 "Application 
of Satellite Radiation Budget Measurements to Study Climate Processes 
of the Earth". Under this cooperative agreement Mr. Shi-Keng Yang, a 
graduate student in Atmospheric Science of the University of Michigan, 
worked with Dr. 6. Louis Smith of the Atmospheric Sciences Division 
at NASA Langley Research Center. 

The research accomplished consists of two parts: 

(a) The development of a model which used climatological data 
for the Earth-Atmosphere system for the calculation of earth 
emitted radiation. This was the Ph.D. dissertation of Mr. 
Shi-Keng Yang, "An Earth Outgoing Longwave Radiation Climate 
Model", which constitutes the main part of this report. 

(b) A study of atmospheric temperature lapse rate. The result 
of this study are summarized in a paper "Further Study on 
Atmospheric Lapse Rate Regimes", which has been submitted 
to the Journal of Atmospheric Sciences. An initial version 
of this paper appears in the semi-annual report on this 
project which was submitted in August, 1984. 


Fred L. Bartman 
Project Director 


ABSTRACT 


AN EARTH OUTGOING LONGWAVE RAO I ATI ON CLIMATE MODEL 

by 

Shi-Keng Yang 

Chairmen: Fred L. Bartman, G. Louis Smitn 

An Earth outgoing longwave radiation (OLWR) climate model has been 
constructed for radiation budget study. The model consists of the 
upward radiative transfer parameterization of Thompson and Warren 
(1982), the cloud cover model of Sherr et (1968) and a monthly 
average climatology defined by the data from Crutcher and Meserve 
(1971) and Taljaard et aj_. (1S69). Additional required information 
is provided by the empirical lOOmb water vapor mixing ratio equation 
of Harries (1976), and the mixing ratio interpolation scheme of Brieg- 
leb and Ramanathan (1982). Cloud top temperature is adjusted so that 
the the calculation would agree with NOAA scanning radiometer measure- 
ments. Both clear sky and cloudy sky cases are calculated and dis- 
cussed for global average, zonal average and world-wide distributed 
cases. The results agree well with the satellite observations. 

The clear sky case shows that the OLWR field is highly modulated 
by water vapor, especially in the tropics. The strongest longitudi- 
nal variation occurs in the tropics. This variation can be mostly 
explained by the strong water vapor gradient. Although in the zonal 
average case the tropics have a minimum in OLWR, the minimum is essen- 
tially contributed by a few very low flux regions, such as the Amazon, 


Indonesia and the Congo. There are regions in the tropics such that 
their OLWR Is as large as that of the subtropics. In the high la- 
titudes, where cold air contains less water vapor, OLWR is basically 
modulated by the surface temperature. Thus, the topographical heat 
capacity becomes a dominant factor in determining the distribution. 

Clouds enhance water vapor modulation of OLWR. Tropical clouds 
have the coldest cloud top temperatures. This again increases the 
longitudinal variation in the region. However, in the polar region, 
where temperature inversion is prominent, cloud top temperature is 
warmer than the surface. Hence, cloud has the effect of increasing 
OLWR. The implication of this cloud mechanism is that the latitudinal 
gradient of net radiation is thus further increased, and the forcing 
of the general atmospheric circulation is substantially different due 
to the increased additional available energy. 

The analysis of the results also suggests that to improve the 
performance of the Budyko-Sel lers type energy balance climate model 
in the tropical region, the parameterization of the longwave cooling 
should include a water vapor absorbing term. 
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CHAPTER 1 


INTRODUCTION 


1.1 Background 

The Earth radiation budget characterizes climate. Net radiation, 
the difference between the absorbed solar radiation and earth emitted 
outgoing longwave (LW) radiation, defines the regions of energy 
sources and energy sinks. These are the dynamic force driving the 
atmosphere-ocean system. An understanding of the radiative processes 
and distribution of the radiation energy is important to gaining 
insight into weather and climate processes, and hence is important 
to improving our prognostic ability, which is one of the most funda- 
mental problems in atmospheric science. 

Earth outgoing LW radiation is an indication of Earth-atmosphere 
conditions. Its temporal and spatial variations result from changes 
in the Earth-atmosphere properties. For example, Lau and Chan (1983) 
used LW as an indicator for studying teleconnections; Gill and Ras- 
musson (1983) and Quiroz (1983) used the LW anomaly to analyze the 
latest " El Nino'* event. 

The introduction of the time-space dependent radiation field into 
numerical prediction schemes will hopefully increase predictability in 
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extended weather forecasts (Wang et al_. , 1983). The study by Belov 
et aK (1971) incorporating radiative heat influx in numerical 
forecasting is another example. In the case of short-term climate 
prediction, radiation maps provide one more set of sample data of 
present day statistics for anomaly analysis, which appears to be 
the only method used in current operational climate forecasts. Many 
theoretical energy balance climate models are simplified to zonally- 
averaged 1-0 models for econony of computation. Although this 
energy balance concept has been widely used in theoretical climate 
models, e.g. Budyko (1969) and Sellers (1973), it has not been 
applied to real time forecasts except for the work of Adem (1969 and 
1981), who used an energy balance model to forecast sea temperature 
anomalies and drought with limited successes. This simplification 
filters out longitudinal variations, including the so called monsoonal 
circulation, hence some interesting pheonomena of the outgoing LW 
radiation features vanish. On the other hand, the top of the hierarchy 
of climate models, GCMs, are more dynamically oriented, but the 
radiation budget has not been well simulated yet (Herman, 1981). 

Radiation budget measurements have long been performed at ground 
bases. London (1957) summarized zonal ly averaged radiation components 
in this manner. It was not until the Explorer 7 satellite (Suomi, et 
al . , 1960) that it was possible to directly measure the earth-space 
radiation exchange, regardless of its crude instrument resolution. 
The Nimbus 6 Earth Radiation budget (ERB) Experiment started a con- 
tinuous time series of measurements in July 1975, which were later 
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succeeded by Nimbus 7 ERB. Bess, et aU (1981), using the method 
developed by Smith and Green (1981), analyzed the LW from the first 
year of Nimbus 6 ERB data thereby revealing some interesting intra- 
annual variations. This method enhances the instrument resolution 
of the wide-field-of-view radiometer. Campbell and Vonder Haar (1980) 
averaged data from some satellite observations and formed a climato- 
logy of radiation budget measurements, which provides very broad in- 
formation on all the earth radiation components. Some of the radia- 
tion features described by these authors have been qualitatively 
explained. The questions on how the radiation field affects weather 
and climate still need to be explored. Beginning in October 1984, 
the Earth radiation budget experiment (ERBE) onboard the ERBE satel- 
lite and NOAA-F satellite will continue to update the measurements 
with more advanced instruments and data analysis procedures. 

Bartman (1980 and 1981) first built a radiation budget oriented 
global albedo climate model. This 2-0 model empirically describes 
the solar radiation transfer process on the earth surface and in the 
atmosphere and calculates planetary albedo at the top of the atmo- 
sphere for each 10° x 10° grid interval. This model can also 
serve as a reference for the climate model in specifying the absorbed 
solar radiation of each region, the radiative heating. No similar 
work, however, has been done on earth emitting LW radiation, the 
radiative cooling. An Earth outgoing longwave radiation climate 
model of this kind will be important to the understanding of our 
Earth-atmospnere climate. 
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1.2 Current Work 

To simulate the earth's outgoing LW radiation field, a 2-D longi- 
tude-latitude varying outgoing LW radiation climate model Is constru- 
cted. The seasonal variation of LW radiation due to variation In 
surface temperature, cloud cover, cloud top temperature and moisture, 
the most Influential elements, will be studied. The model also 
provides a tool to study some Interesting features of the LW field, 
such as the LW maximum in the subtropics and drastic longitudinal 
variations along the tropics. 

The model consists of a radiation calculation and a data base of 
climatologies of the input parameters, i.e. the monthly average fields 
of surface temperature, cloud cover, cloud top temperature and mois- 
ture over the globe. For efficiency, a parameterization by Thompson 
and Warren (1982), which is derived from Wiscombe's radiative transfer 
model (1976), is used for the radiation calculation. This will be 
discussed in Chapter 2. 

For describing climatologies of input parameters, a 10° x 10° grid 
system is selected. A climatological data set from Crutcher and 
Meserve (1970) and Taljaard et al_. (1971) is used to provide temper- 
atures, dew point and geopotential height for various altitudes. 
This will be discussed in Chapter 3. 

In Chapter 4, because of the lack of data at 300mb and 200mb, 
the interpolation method used by Brlegleb and Ramanathan (1982) is 
incorporated to get 300mb and 200mb specific humidities with the 
formula derived by Harries (1976) for lOOmb. Also a theoretical 
limatology of clear sky outgoing LW radiation is presented. 


In Chapter 5, Sherr et al_. (1968) cloud type classification and 
cloud cover data are incorporated with satellite observations to 
derive a scheme for determining cloud top temperature. In Chapter 
6, the calculated cloudy sky outgoing longwave radiation is discussed 
in detail. Finally, the conclusions and recommendations are stated 
in Chapter 7. In most cases, only the results of the extreme seasons, 
January and July, are discussed. The cases of the transient seasons, 
April and October, are generally more moderate for all the features. 
These results are organized in the Appendices for comparison. 

Fig.l summarizes the methodology of this outgoing LW radiation 
climate model . This model explains quantitatively each parameter's 
contribution to the LW field, and thus gives a better understanding 
of the physics of longwave radiation in the earth climate. 
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FIGURE I . METHODOLOGY OF EARTH OUTGOING LONG WAVE Cl MATE MODEL 
RH a bright vrightrd mr«n rrUtivr humidity , m chrptrr 2. 
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PARAMETERIZATION OF EARTH OUTGOING LONGWAVE RADIATION 

2.1 Introduction 

As briefly mentioned in the previous chapter, an Earth outgoing 
longwave radiation (OLWR) parameterization by Thompson and Warren 
(1982, hereafter denoted as TW) is used for this study. In section 
2.1, the reason for choosing this parameterization is discussed. 
In section 2.2, the different parameterizations of OLWR are reviewed, 
and in section 2.3, the TW parameterization is outlined and discussed. 

2.1.1 Longwave Parameterization and Climate Model 

In spite of the fact that the basic radiative transfer equation 
is simple for the plane-parallel atmosphere, the solutions are not 
straightforward for the reasons that (1) the atmosphere is inhomoge- 
neous and (2) the molecular absorption spectra are extremely compli- 
cated. Often a band model is used to alleviate the cumbersome line- 
by-line integration of absorption for transmission functions. Also 
the Curtis-Godson approximation is widely used to replace an inho- 
mogeneous optical depth with an equivalent homogeneous one. Never- 
theless the selection of bandwidth and the treatment of the distri- 
bution of line strengths are somewhat arbitrary, e.g. band width 
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could be as narrow as lOcnr 1 for a narrow band model to more than 
100cm" 1 for a wide band model. The classic El sasser regular band 
model and Goody random model are the backbone of most derivations. 

To verify band models, a detailed line-by-line calculation Is 
Indispensable. Drayson (1967) first performed the calculation with 
very fine line width of 0.001 cm" 1 over line center and 0.1 cm* 1 
over line wings and integrated over very thin layers so that the 
inhomogeneous path could be accounted for. 

McClatchey et al_. (1972) compiled absorption line dat*., from 
microwave to the visible, not only for COg, O3 and H2O bands but 
as well for the other trace gases which are radiatlvely active In 
the earth atmosphere. This set of absorption line data is being 
updated constantly and currently contains 159,000 lines. It is the 
most complete data tape available. McClatchey et aj_. also developed 
a line-by-line calculation program for transmisson and radiation, 
FASC00E. The main problem In intercomparing detailed line-by-line 
models with band models using published results is that different 
models have different Input conditions and different resolutions 
(Luther, 1983), which makes direct comparison with output results 
meaningless. Luther has conducted a workshop in 1984 to compare 
the results from different models. Some of the preliminary results 
have shown that great uncertainty on the water vapor continuum absorp- 
tion could cause significant error. This work will substantially 
improve our understanding of the performance of radiation models. 

A detailed radiative transfer model is expensive to use. Even a 
one-dimension Quasi-Random model takes a few seconds computation 
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time. To calculate a 10 0 x 10 0 grid global distribution 648 grid 
points Is very costly. A well parameterized radiative scheme can 
give results as good as the radiative transfer model. Also, para- 
meterization creates an Indispensable Interface with available input 
data. 

Based on these arguments, a parameterization is selected for this 
study to calculate outgoing LW radiation directly. 

2.2 A Review of Earth Outgoing Longwave Radiation Parameterization 
Budyko (1969) empirically parameterized the Earth's outgoing LW 
radiation using monthly data of 260 stations as 

F- a+bT-(ai+biT)C 2.1 

Where F is outgoing LW radiation at the top of atmosphere in w/m^ 
T is the surface temperature. 

C Is cloud cover. 

a,b,ai and b^ are constants. 

This formulation was compared with the result of the theoretical 
radiative transfer model of Manabe and Wetherald (1967) and claimed to 
be in good agreement for the clear sky case but with some differences 
for the cloudy sky case (Rudyko, 1969). 

Subsequently, a series of experiments was performed on the 
derivation of coefficients using either an empirical approach or a 
theoretical model calculations, such as Schneider (1972), Cess (1974) 
and Gupta et a_k (1978). A summary of these coefficients are tabulated 
in Table A. 18. It is clear that they are in good agreement for the 
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global average case. 

As pointed out by Thompson and Warren (1982), the disadvantages 
of empirically derived parameterzations are: 

(1) the coefficients vary from data set to data set. Currently 
there is no generally recognized best value of b for equ.2.1. 
Recently Short et al_. (1984) have shown that b is temporally and 
spacially dependent in the tropics. 

(2) When two physically distinguishable but highly correlated 
factors, such as cloud and humidity, both influence F, it is very 
difficult to separate them and establish their relative magnitude 
by regression analysis. 

The advantage of the empirical approach is that the parameterization 
might implicitly account for some unsuspected feedback effect in the 
climate system that might otherwise be neglected in a climate model. 
On the other hand, the advantage of the parameterization derived 
from the transfer model is that the relation between input atmos- 
pheric condition and the corresponding output F are perfectly known. 
The disadvantage is that the result would possess the defects of 
the model. For example, the concentration of stratospheric water 
vapor is not well known yet (see Chapter 4). In most radiative 
transfer models it is assumed to be about 3ppm. Altering strato- 
spheric water vapor concentration would change the longwave opacity 
and thus the outgoing LW flux. Wang et aU (1976) estimated that 
doubling water vapor concentration would decrease LW flux by about 
2W/m _ 2. Therefore, a model could have a potential error due to 
the difference between the assumption and the true value. A parame- 
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terlzatlon derived from such a model would thus Inherit this error 
as wel 1 . 

Different formulations of the outgoing LW radiation have been 
tried by other investigators, especially to avoid the uncertainty of 
cloud amount. Sellers (1969) used an atmospheric attenuation coeffi- 
cient to replace cloudiness in a simple zonally averaged climate 
model. Fromm (1982) and Jensenius ^et ^al_. (1978) derived a parame- 
terization using other meteorological parameters, such as pressure 
tendency, wind magnitude, vortlcity, etc., which cloudiness is impli- 
citly related to. This approach is relevent for short term numerical 
model applications. 

Since the intriguing question of how the temporal and spacial 
variation of water vapor Impact on the outgoing LW radiation is the 
focus of this study, the Thompson and Warren (1982) parameterization 
is chosen for its special feature of water vapor as an input predictor. 
This allows one to study the relative impacts on the outgoing LW 
radiation of both cloud and water vapor. For example, in the zonally 
averaged case, the outgoing LW radiation exhibits a minimum near the 
equator 20 to 40 W/m^ smaller than the peaks at the subtropics. In 
a calculation by Warren and Thompson (1983), they found that one 
third of the dip can be explained by the variation of atmospheric 
water vapor. Only a small fraction is attributed to cloudiness, if 
the published cloud climatology is used. The rest remains unexplain- 
ed. They thus suspect that the cirrus and cumulonimbus amounts near 
equator are much higher than the published cloud climatology. In 
chapter 4, we will further discuss how the atmospheric water vapor 
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modulates outgoing LW radiation. 

A brief discription of the important features of this paramete- 
rization follows, as well as some crucial information not mentioned 
in the original work. 

2.3 Thompson and Warren Parameterization 
2.3.1 Wiscombe Model 

In Thompson and Warren (1982), outgoing LW parameterization is 
derived from a radiative-transfer model. The model used is Wiscombe 
(1975) ATRAD, which employs the exponential -sum-fitting of transmi- 
sion function method. Essentially the transmission functions are 
those of LOWTRAN 5 (Kneizys et al_. , 1980). The LW spectrum is covered 
from 3 to 500 y m, but differs from LOWTRAN 5 in that water vapor 
continuum is included. This is important in the tropics with its 
high water vapor content; omission of the continuum absorption in 
the radiative trailer calculation in early GCMs led to the calcula- 
tion of very low atmospheric temperatures. 

2.3.2 Assumptions 
1) Major Gases: 

Mixing ratios of the gases, COg, N 2 O, CH 4 , CO, N 2 and O 2 remain 
constant at all altitudes at the following values: 320, 0.28, 1.6, 
0.075, 7.965xl0 5 and 2.095x10^ ppm respectively. Current prominent 
attention on the climatic impact of anthropogenic CO 2 increase 
leads to the question of how outgoing LW radiation responds. Increa- 
ing the concentration of CO 2 in the atmosphere will result in: 
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(1) an increase of tropospheric-surface temperature because of the 
greenhouse effect. This can be measured from the net change in LW 
flux at the tropopause (Ramanathan et al_. , 1979). 

(2) Cool stratospheric temperature resulting from increasing strato- 
spheric opacity, i.e. more effective emission. 

However, the influence on the outgoing LW at the top of the atmo- 
sphere is roughly a factor of 2 smaller than the net change of LW at 
the tropopause. Doubling the CO2 concentration will only change 
outgoing LW radiation at the T.O.A. about 2 w/m 2 . Error from 
neglecting the variation of CO2 is smaller than the errors from 
other sources, such as data error and cloud fraction. A broad band 
channel which measures the total outgoing LW is not sensitive enough 
to detect such a small fractional variation. However an instrument 
on a satellite using a narrow channel centered at the CO2 absorption 
band can detect more significant change in that spectral range due 
to atmospheric CO2 change {Kiehl, 1983 and Charlock, 1984b). 

2) Ozone: 

Ozone O3 profiles are interpolated from the McClatchey (1972) 
model atmosphere. In the LW spectrum, O3 is only significant at the 
9.6 jum band. Changing the O3 concentration causes significant 
change in the outgoing flux of this band only; for broad band the 
total change is very small. For example, a change of O3 concentra- 
tion from the tropical atmospheric value (.246 atm cm) to its mid 
latitude atmospheric value (.345 atm cm), changes the outgoing flux 
for the 9.6um band (1000 to 1120cm“l) about 3.7w/m 2 , but the 
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whole LW spectrum (400 to 1400cm* 1 ) only changes about 4.3 w/m 2 
(Charlock, 1984a). In the cloudy sky case, fixing the O 3 distribu- 
tion should have less error. 

3.) Aerosols: 

The question of the possible impact of increasing aerosols on 
climate is gaining as much attention as the CO 2 greenhouse effect 
these days. Some inconclusive conjectures predict that anthropogenic 
aerosols' impact will be of the same order as that of CO 2 in the 
future. Nevertheless, modeling aerosols is far more complicated 
than CO 2 (Ramanathan and Coakley, 1978). Unlike CO 2 , aerosols are 
a more localized problem, which do not diffuse homogeneously to 
the whole atmosphere. In the United States, current aerosol concen- 
tration in urban areas can be three times greater than those over 
rural areas. Fallout time of normal tropospheric aerosols is 1 to 2 
weeks, while that of stratospheri c aerosols can be as long as years. 
Thus the concentrations, size distributions (.01 to 20 um) and opti- 
cal properties have high dependence on both time and space. 

Aerosols not only reduce outgoing longwave radiation as Joes CO 2 
or any other atmospheric constituents, but also scatter and absorb 
incoming solar radiation covering the entire spectrum. The scatter- 
ing-absorption processes can result in either increase or decrease 
of planetary albedo depending on the nature of a particular type of 
aerosol. Some absorb (scatter) more solar energy than they scatter 
(absorb), which will decrease (increase) planetary albedo when the 
concentration increases. Hence, the s'ign of heating on the earth 
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surface temperature from Increasing aerosol would be from the diff- 
erence of the two processes: the net planetary albedo change and 
the decrease of outgoing longwave radiation. Coakley and Grams 
(1976) found that both small particles (radii smaller than 0.05 urn) 
and large particles (radii larger than 1.0 -urn) generally have a 
greater influence on longwave radiation than on incoming solar radi- 
ation, therefore these particles contribute to warming the surface 
temperature. Particles of the intermediate size affect the incoming 
solar radiation more strongly than they affect the longwave radiation 
and thus contribute to cooling the surface temperature. Reck (1974) 
and Wang and Domoto (1974) also found that the surface albedo has a 
crucial role in determining the sign of heating. When the surface 
albedo is higher than a critical value, increasing aerosols will 
cause surface temperature warming. However, not all of the study 
results agree. Most of these experiments were carried out with the 
one dimensional radiative convective models or energy balance models 
with the size distribution of spherical particle between 0.1 *jm and 
1 wm, which is not representative in simulating a localized problem. 
One can thus only expect to gain the qualitative understanding of 
aerosol physics. 

In conclusion, these studies showed that the aerosol self compe- 
ting effects might have important climatic influence on the surface 
temperature. For the purpose of the present investigation on outgoing 
LW radiation, the effect of aerosols will be neglected. 
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4) Surface emissivity: 

Surface emissivity is assumed fixed at .95 by the TW parameteri- 
zation. Data of global distribution of surface emissivity are not 
available yet. In the 8 to 14 urn region, which accounts for most of 
the surface radiation contributiou to outgoing LW radiation, this is 
the intermediate between extremes of .90 for sand desert and .99 
for snow. Ocean is about .98. Outside this window band the varia- 
tion of surface emissivity has little effect on the outgoing LW 
radiation from the surface. TW showed that a change of .05 in 

surface emissivity will vary clear sky outgoing LW radiation by 
only 1 to 3 W/m^ depending on atmospheric opacity. 

5) Discontinuity of temperature at the surface: 

Surface and surface-air temperatures are assumed equal. In fact 
there is usually a few degrees Kelvin difference of diurnal mean 
temperatures between surface air and ground temperature. Tests show 
that a one degree temperature difference between surface-air and 
ground would contribute a .7 to 1.5 W/m^ change in clear sky out- 
going LW (Thompson and Warren, 1982), which means that with thi c 
assumption the model can not simulate well the LW diurnal variation. 
For monthly averaged cases, this assumption has little effect on 
monthly mean outgoing LW. 

6) Surface Elevation. 

Surface elevation is not included as a predictor of F. In the 
case of low temperature, such as in the Antarctic, this assumption 



contributes little error because there is less water in the air 
column. However, if the surface temperature is higher, such as in 
Tibet, the error can be higher. TW tested this by changing surface 
pressure from 1013mb to 700 mb, which can increase outgoing LW by 
8w/m 2 if the 40°N zonal average temperature is used. However, if 
the local temperature, which is lower than zonal average temperature, 
is used, the difference will substantially decrease. For very high 
elevation, the surface temperature is generally lower, which will 
reduce water vapor content, and thus produce less error. 

2.3.3 Predictors and Characteristics of Parameterization 
The parameterization now has the form: 

F=C 1 -C 2 A c 2.2 

where and C 2 are coefficients, and are now functions of surface 
temperature, cloud top temperature and height weighted mean relative 
humidity which will be further discussed later. A c is effective 
cloud fraction, Cj accounts for the clear sky ougoing LW, and C 2 A c 
is the modification term which accounts for the effect of cloudiness 
on the outgoing LW. 

A more physically straightforward form of equ 2.2 is 

F-Ci(1-A c )+C' 2 A c 2.3 

where C‘ 2 *Ci-C 2 . This equation simply separates the clear sky part 
from the overcast part. However, equation 2.2 is more often used 
because C 2 = dF/dA c measures the cloud LW feedback (Cess and Ramana- 
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than, 1978). 

For outgoing LW budget on a global scale, the surface radiation 
contributes about 10 to 15%, the atmosphere about 24% to 54% and 
the remaining is cloud radiation (Wallace and Hobbs, 1977). The 
infrared cooling of the troposphere is highly dependent on the temper- 
ature at different levels of the atmosphere. Since the temperature 
sounding profiles at most latitudes have similar characterise cs 
(Liou, 1980), using surface temperature to represent both surface 
and atmospheric LW emission is a very relevant and convenient simpli- 
fication. As altitude increases, the absorbent contents decrease, 
and the atmosphere above becomes more transparent. Above the cloud 
top, generally it is clear enough for cloud radiation to escape to 
outer space very efficiently. Thus, using cloud top temperature as 
another predictor is well justified. 

Another predictor, effective cloudiness, is a simplification of 
compressing multilayer clouds into one simple black cloud layer. In 
reality, this assumption may be valid for single layer low thick 
clouds. The high thin cirrus type clouds generally exhibit spectral 
emissi vities smaller than .5. In those cases without overlapping, 
the effective cloudiness would be emissivity times cloud fraction 
A c , such that a grid area is divided into a simple clear part and 
an overcast part. This is a very idealistic assumption. However, 
until one knows more about complex cloud behavior, this is a feasible 
approach in terms of modeling. Clouds will be further discussed in 
Chapter 5. 

Atmospheric water vapor is the most efficient LW absorber. As 
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observed from the space, one would find that atmospheric water vapor 
is the only LW sink compared to the other three predictors, T s , T c 
and A c , just discussed above. In the earlier empirical LW parame- 
terization, this term is often omitted probably due to lack of simul- 
taneous observation. As will be discussed in Chapter 4, even in 
this modern age there are still some substantial problems in the 
measurement of atmospheric water vapor. Yet there are many ways to 
describe atmospheric water vapor, namely relative humidity, dew 
point, mixing ratio, specific humidity or precipitable water. For 
the TW parameterization, relative humidity is chosen from a pertur- 
bation test on the McClatchey atmospheres (McClatchey et al^, 1972). 
They changed the tropospheric temperature lapse rate, which is not 
included as a predictor, in Wiscombe's radative transfer model, 
and used different water vapor parameters to calculate outgoing LW 
flux. They found that the results are not sensitive to lapse rate 
change if relative humidity is used. Therefore they defined a height 
weighted mean relative humidity as: 


RH = 


TTki 



RH(z) dz 


2.4 


This predictor gives greater weight to relative humidity at high 
altitudes, because high altitudes are more transparent to space than 
the lower levels and therefore contribute more to F. Using 12km as 
the upper limit for integration is mainly due to the convenience of 
parameterzation. The altitude of 12km is approximately the average 
height of the tropopause. This simplification eliminates the latitu- 



dinal and seasonal variation of the tropopause height, which is 
above 12km at the tropics and below 12km at higher latitudes. The 
error of underestimating tropopause height at the tropics and over- 
estimating in the higher latitude Is difficult to measure due to lack 
of observation in this region. This will be further discussed in 
Chapter 4. 

Tables 2.1 and 2.2 from TW show the parameterization of Cj and 
C 2 respectively. Corresponding figures 2.1 and 2.2 illustrate Cj 
as a function of surface temperature and height weighted mean relative 
humidity RH, and C 2 as a function of Rh[ and the temperature difference 
between the surface and cloud top. 


Table 2.1 Parameterization for dear-sky outgoing IR irradiance 
at the top of the atmosphere (z « 50 km) after TW. 
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Table 2.2 Parameterization of cloud modification term for out- 
qolnq IR irradlance at the top of the atmosphere (z - 50 km) after TW 


c 2 A c • cloud modification term (W m* 2 ) [see Eq. (2.2)]. 
Ac ■ fractional total cloud amount. 

T$c ■ T s * ^c < ^sc * 60K). 

T s « surface air temperature (-118°C < T c < 57°C). 

T c - cloud top temperature (-118°C < T c < 57°C). 

RH « height-mean relative humidity [Eq. (2.4)] 

(0.2 < RH < 1.0). 

PARAMETERIZATION: 

C2 * C 1 (T s . RH) ” c l( T C» RH) + f(T sc» M) 

(for cj see Table TTl) 

f(T sc .M) * f l(Tsc) + Wso™) 
fj * d 0 + d 1 T sc + d 2 T 2 sc 

f 2 * d 3 T sc( T sc + d 4) (M + d 5) 

do * -3.1 

dj * -0.4146 

d 2 » 4.084 x 10" 3 

d 3 - -4.44 x 10* 3 

d4 ■ 80.0 


d 


-0.40 
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Clear 

sky 

outgoing 
IR at T.O.A. 
(W/m- 2 ) 



Surface temperature (°C) 


Figure 2.1 


Family of curves depicting the parameterization of clear 
sky outgoing longwave irradlance las a function of T s and 


RH. 
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Parameterization of C 2 

RH = 30% Radiative 
60% + Transfer 



Figure 2.2 The cloud modification parameter C? as a function of 
T $ -T c and RH for McClatchey et al_. (1972) tropical 
temperature profile. 
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CHAPTER 3 


DATA AND DATA MANAGEMENT 

A climatological monthly mean data tape supplied by NCAR provides 
a great portion of the input data base for this study. In sec.3.1, 
we briefly discuss the content of the tape. In sec. 3. 2, the data 
averaging is mentioned. 

3.1 Data 

The data tape consists of two data files, one for each hemisphere. 
Both files have data on a 5°x5° horizontal grid system. Meteorological 
variables on the tape are: temperature, dew point, geopotential 
height, standard deviation of these parameters and geostrophic wind. 
On the surface, sea surface pressure replaces geopotential height. 
Dew points are only recorded from the surface through 500mb. Two 
different forms of data are provided, one is hand analyzed, and 
the other is smoothed data. The smoothed temperature, dew point and 
geopotential height are selected for this study for the reason that 
the climatic fields are smooth after long term averaging. 

Crutcher and Meserve (1970) used the data gathered by the World 
Meteorological Organization (WMO) climate data program for the North- 
ern Hemisphere. The period of surface data is from 1931 through 1960, 
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and the period of upper air data covers the years from 1950 through 
1964. However, many stations have much less data; some only have 
data 2 to 3 years. The data for the Southern Hemisphere were compiled 
by Taljaard, et jj 1_. (1969). Data sources are mostly from the WMO 
International Geophysical Year, Air Force Air Weather Service and 
historical data from various sources. Station density is less in 
the Southern Hemisphere due to more ocean coverage. Isopleths were 
subjectively drawn across void grid points, and therefore the relia- 
bility of each grid varies depending on the availability of data. 

Data are also subject to instrumentation errors, which include 
bad calibration and different performance from different instrument 
manufacturers. In the troposphere, the standard deviation of tempera- 
ture measurement from using different instruments is about .5°C, and 
in the stratosphere is about 2°C. Errors are also introduced in 
communication procesres, such as radio noise, teletype errors and 
human error. 

Computer routines were designed by Crutcher and Meserve (1970) 
and Jenne, et aK (1974) to eliminate the unreasonable data. All of 
the errors mentioned above are summarized in the variance analysis 
and tabulated ?s standard deviations for each parameter except geos- 
trophic wind, which we did not use. 

3 2 Data Averaging 

A simple scheme is defined to merge the two files of the Northern 
Hemisphere and the Southern Hemisphere together, and to average the 
orginal 5°x5° grid system to 10°xl0°. The reasons for using 10°xl0° 
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grid system is that the cloud data are not available in more detail. 
Also, the Thompson and Warren (1982) parameterization is limited to 
coarse resolution application. 

The new grid points are selected at latitude 85°N, 75°N,..., 75°S, 
85°S on longitude 5°W, 15°W, ..., 355°W. Data from all nine grid 
points around each new grid point (see Fig 3.1) are included in the 
averaging process. However, near the poles, i.e. 85°N and 85°S, 
the original data at 90° are only counted once to avoid overweighting 
the repetition of same data at every longitude. This procedure 
also smooths the data slightly. 



Figure 3.1 Grid points for the new 10°xl0° system relative to the 
NCAR 5°x5° system; (a) for ordinary latitude, (b) at 
the poles. 



CHAPTER 4 


WATER VAPOR ANO CLEAR SKY OUTGOING LONGWAVE RADIATION 
4.1 Introduction 

As mentioned in Chapter 2, one of the parameters used in the 
longwave outgoing flux calculation is the average relative humidity, 
defined as: 


, J2km 

Bit* TFkm ( RH ( 2 > dz 4 - 1 

where the pressure at 12km height is about 190mb (U.S. Standard 
Atmosphere, 1976). To implement equation 4.1, one needs (1) suffi- 
cient water vapor data from the surface to 12km, (2) an algorithm 
to convert dew point or mixing ratio to relative humidity, and (3) 
an integration scheme. However, the observational data set described 
in Chapter 3 only provides dew point information from the surface to 
500mb. This data between 500mb and 12km must be furnished through 
another source. To fulfill this, we adapt an empirical formula by 
Harries (1976), which describes the water vapor distribution at lOOmb, 
and an interpolating scheme which incorporates the Harries empirical 
formula and the observational data to generate water vapor data at 
200mb and 300mb. 
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In section 4.2, we will d’^cusi the atmospheric water vapor. In 
section 4.3, the interpolation scheme will be discussed. In section 
4.4, the method of converting dew point temperature and mixing ratio 
to relative humidity Is discussed. Section 4.5 describes the integra- 
tion of equation 4.1. In section 4.6, the world wide RH and clear 
sky outgoing longwave radi;tion are presented. In section 4.7, 
the sensitivities of th" interpolation scheme and averaging scheme 
are discussed. 

4.2 Atmospheric Wate r Vapor 

Although it 4 s , jt necessarily true that atmospheric water vapor, 
when expressed as relative humidity, decreases with altitude, it is 
well recognized that mixing ratio decreases rapidly with altitude 
(Smith, 1966) as it is an exponential function of temperature. Unable 
to cope with such a drastic change, the hygrometers used in conven- 
tional radiosondes generally lose their sensitivities as temperature 
goes below -40°C. Even under the warmest atmospheric conditions, 
radiosonde humidity observations seldom exceed an altitude of 10 km 
(U.S. Standard Atmosphere, 1976). London (1957) urged water vapor 
measurement in the upper troposphere for its importance to the radi- 
ation budget. Yet, Routhier and Davis (1980) analyzed the data 
collected from the Global Atmospheric Measurements Experiment of 
Tropospheric Aerosols and Gases, and suggested that some of the 
global climatological dew point data in the middle troposphere 
measured from the conventional method was too high due to such an 


instrumental limitation. 
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Meanwhile, several in situ stratospheric water vapor experiments 
with special instruments, which included part of the higher tropo- 
sphere had been performed. The most well known are the long term 
series of the measurements conducted by Mastenbrook (1968). These 
measurements used balloon soundings with a frost-point hygrometer. 
The results showed that there is a seasonal variation in mixing 
ratio with the lowest mixing ratio in late winter and early spring 
and the highest in late summer and early fall. Also, the magnitude 
of the seasonal change decreases wit. height. He thus suggested 
that the temperature of the tropopause regulates the level of strato- 
spheric moisture. 

However, not all of the experiments agree well with each other. 
Figure 4.1, from a review of Harries (1976), shows results from the 
major measurements performed at different times, seasons, latitudes 
and with different instruments. By rejecting some of the conflicting 
data, he managed to formulate an equation governing the latitudinal 
variation at lOOmb, but without the time variation: 

W < J > * 3 * i°' 6 /n + ) 4.2 

where J is the index for latitude; 1, 2, 3, ...» 35 correspoing to 
-85°, -80°, ..., 85°. This equation was later implemented by 
Briegleb and Ramanathan (1982) for their clear sky albedo calculation. 

A more recent review on stratospheric water vapor was done by 
Ell sasser (1983). He tried to interpret the tropospheric-strato- 
spheric water vapor transfer mechanism with minimum rejection of 
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Figure 4.1 Summary of observations of vertical profiles of high 
altitude water vapor compiled by Harries (1976). 
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conflicting observational data. Basically* the stratospheric water 
vapor is flushed by air entering only through the tropical tropopause 
and leaving only through the polar tropopauses or tropopause gaps. 
There is also a possible chemical source, methane oxidation, which 
however, does not have significant effect. 

In contrast to the in situ measurements, the recent LIMS Experi- 
ment on the Nimbus 7 satellite measured the water vapor distribution 
higher than lOOmb. This experiment revealed very remarkable strato- 
spheric water vapor information (Russell et al_. , 1984): there is a 
broad tropical minimum in mixing ratio, extending from lOOmb to 
about 30mb as shown in Figure 4.2, which coincides with the "cold 
trap", the coldest temperature region in the stratosphere. The data 
also showed a poleward gradient at all levels up to 4mb. These 
results are very different from the perceptions of the previous 
investigations and might have very substantial implications for 
water vapor budget and stratospheric general circulation studies. 
In the process of the present study, this result was not available. 

More than 25 years after London (1957), the climatological water 
vapor distribution between 500mb and lOOmb over the globe is still 
unavailable except for some zonally averaged data for the Northern 
Hemisphere from London (1957) and from Oort and Rasmusson (1971). 
The latter raises some problem in converting mixing ratio to relative 
humidity, which will be discussed in section 4.4. Some data for the 
tropics from Newell, et a]_. (1972) are also available. Clearly, 
there is a demand for refining current sounding instruments in order 
that the measurement be operable under low water vapor pressure, low 
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pressure conditions, so that compiled climatologies such as Thompson 
and Warren’s (1982) zonal 1y averaged relative humidity or the results 
from complicated general circulation models, such as Sirutis, et al. 
(1980) can be verified. 

4.3 Vertical Interpolation for 300mb and 200mb Specific Humidity 
from 500mb and lQOmb 

As mentioned in the section 4.1, we need to fill in water vapor 
data higher than 500mb in order to perform the vertical integration. 

An interpolation method similar to that used by Briegleb and 
Ramanathan (1982) is adapted for this purpose. This method incorpo- 
rates Harries (1976) empirical formula for lOOmb water vapor, equation 
4.2, and the NCAR data base described in Chapter 3 as boundary condi- 
tions. It will be shown in section 4.7 that the error of the final 
longwave flux calculation due to this interpolation is minimal. 

Briegleb and Ramanathan studied the radiosonde measurements from 
Rogers and Walshaw (1966) and found that the variation of the logarithm 
specific humidity 1 n (Q ) with pressure P from 500mb to lOOmb has a 
quadratic relationship as 1 n (Q) ~ p2. Thus, if the data of these 
two layers are available, the data of intermediate layers can be 
obtained through interpolation. A special interpolation function 
with a single parameter is selected such that only upper and lower 
boundary conditions are needed, compared to polynomial interpolation 
which generally needs more than three data points specified. The 
interpolation formula is 


in Q(P) - £ (in Q 100mb + In Q 500mb ) + 


(ln Q 500mb ' ,n QlOOmb^^r 1 + 4 * 3 

«(l-n)( 1n Ql00mb \ ^ Q 5°On~b)(P j 300 ) ; 

where subscripts indicate pressure levels and n is an interpolation 
factor. When n=l, this is a linear interpolation and when n is adjus- 
ted to be less than 1, Q can be fitted to the quadratic curves re- 
quired for climatological specific humidity, see Figure 4.3. 

In order that the zonal average Q's of 300mb vary the same way as 
the observations of Oort and Rasmusson (1971), and for the reason that 
Q varies from the tropics to the poles only by a factor of 3 at the 
lOOmb pressure level, while at the 500mb pressure level it varies by 
a factor of 10 (Briegleb and Ramanathan 1982), n must be latitu- 
dinally dependent so that there is more linear interpolation in the 
poles and more quadratic interpolation in the tropics. Both can be 
accomplished by letting n, in equation 4.3, vary as a function of 
latitude as: 


2 

n(j) - 0.86 [1 * ] 4.4 

where j is the index of latitude, 1,2,3,... referring to latitudes 
-87.5°, -82.5°, ..., 87.5°, respectively. 

To restrict overestimation when compared to the climatology of 
Rogers and Walshaw (1966) another constraint needs to be applied: when 
equation 4.4 gives n(j) > .92, it is set to .92. 

Using this method to calculate relative humidity will result in 
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Water vapor mixing ratio 

Figure 4.3 Characteristics of the interpolation function equation 
4.3 for different values of parameter n. 
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a few points with a value higher than 100% at 300mb and 200mb. For 
those cases, the value is set equal to 100%. This Is not unexpected 
for coupling two different data sets together. Even using matched 
zonal average specific humidity and temperature, such as Oort and 
Rasmusson (1971), will yield RH higher than 100%. The reason for 
this will be discussed in the latter part of the next section. 

4.4 Relative Humidity Calculations from Dew Point Temperature and 
Mixing Ratio 

When dew point temperature (T^ ) or mixing ratio (W) are available 
either from measurements or from interpolation as generated from the 
method described in the previous section, one must convert them to 
relative humidity (RH). The definition of RH according to List (1959) 
is 

RH=W/W S 4.5 

where the subscript s denotes saturation and W is defined as 

W*,622 e/(p-e) 4.6 

where e is water vapor pressure. Equation 4.6 can be written 
as 

RH=(e/e s )(p-e s /p-e) 4.7 

The second parentheses on the right hand side approximates to unity, 
so that Thompson and Warren (1982) simplified the definition to 
e/e sw t > 0°C 

RH* 4.8 

e/e si t < 0°C 

where e sw is the saturation water vapor pressure over water, and e s -j 
is the saturation water vapor pressure over ice. There are many 
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methods available for calculation of e sw and e s -|. One can use the 
Integrated form of the Clauslus-Clapeyron equation, the equations by 
List (1959), Richard (1971) or Parish and Putnum (1974). The latter 
three are empirical equations which fit the observed water vapor 
pressure; any one of them can give slightly better results than the 
standard Clauslus-Clapeyron equation In some narrower specific tem- 
perature range. In this study, we follow Thompson and Warren 
(1982) using Richards (1971) for e sw : 

e sw = 1013.25 EXP(13.3l85t-1.9760t 2 -0.6445t 3 -0.1299t 4 ) 4.9 

where t»l-373.15 T _1 and T is the temperature in Kelvin. 

Using List (1959) for e S j : 

L °9l0 e s1 —9.09718(x-l)-3. 56654 Log 10 x+0. 876793 (1-x’ 1 ) 
+Logio6.1071 
where x*273.16 T”l. 

By definition W=e/(p-e), and Q=e/p. Because e is much smaller 
than P, W approximates Q through the meteorological range. In calcu- 
lating relative humidity at the 300mb to lOOmb levels, W is considered 
to be equal to Q, which Is calculated by the Interpolation scheme of 
the section 4.3. Also the corresponding W s must oe calculated and 
transformed to e s . From the surface to the 500mb level, e and e s 
are calculated from dew point and temperature. 

One caution should be taken in selecting the phase change temper- 
ature. Thompson and Warren arbitrarily selected 0°C in their study. 
In the atmosphere, it is very common that water vapor exists in 
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supercooled form. As pointed out by List (1959), water vapor can 
exist in the air temperature below 0°C as supercooled; however, 
hygrometers which essentially only respond to water vapor indicate 
relative humidity with respect to water vapor only, regardless of 
the temperature Thus, -in relative humidity calculations, improper 
selection of the phase change temperature would significantly intro- 
duce error when comparing water vapor pressure to saturation vapor 
pressure of the wrong water phase. To alleviate this error, one 
should select a phase change temperature below 0°C, such as London 
(1957) did. 

Another significant problem in converting climatological specific 
humidity or mixing ratio to relative humidity is that relative humi- 
dity is often overestimated. The error is not directly from the 
converting scheme, but from the time or space averaging scheme for 
the climatology. This can be clearly illustrated by the following 
example in averaging two cases. 

Case I. temperature 273°K, relative humidity 90%, e=5.50mb. 

Case II. temperature 279°K, relative humidity 84%, e=7.85mb. 

So the averaged temperature is 276°K and the averaged e=6.67mb. The 
e $ of the averaged temperature 276°K is 7.56mb. If one uses averaged 
e and the e s based on averaged temperature to calculate relative 
humidity, the result would be 88.2%. However, if one just averages 
the relative humidities, the result would be 87%. A 1.2% difference 
is created. As the sample's deviation increases, the result would 
get worse. This is simply caused by the fact that e s is an exponen- 
tial function of temperature which is concave upward, as shown in 
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Figure 4.4. The water vapor pressure is not supposed to be linearly 
averaged as temperature. Thus when one uses specific humidity and 
temperature for temporal and space averaged climatological data as 
Oort and Rasmusson (1971) and converts them to relative humidity, 
the result is overestimated. In many cases of this data set, the 
results are higher than 100%, (see Table 4.1). On the other hand, if 
one uses averaged dew point temperature to calculate relative humi- 
dity, then it would be underestimated. 

There is no straightforward remedy for this problem unless one 
has sample relative humidity data to start with and directly averages 
over relative humidity. 

4.5 Vertical Integration of relative humidity 

To perform the vertical integration, equation 4.1 is changed to 
numerical form as 

RH= l/12km l (RH i +RH i+1 )(gph i+1 -gph i )/2 4.11 

where gph is geopotential height, subscript i is the index of each 
pressure level from 1 to 6 as surface, 850mb, 700mb, .... 200mb. 

As mentioned earlier, pressure level 200mb is located around 
12km geopotential height, depending on the location. Since equation 
4.1 is an integration function with upper bound at 12km, we rep' ace 
the relative humidity at 200mb by relative humidity at 12km. These 
values are obtained by linear interpolation alor^ the geopotential 
height between 300mb and 200mb or 200mb and lOOmb, depending on 
where the geopotential height of 200mb is located. 
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and Clear Sky Outgoing Longwave Radiation Fields 
4.6.1 Zonal 1y Averaged RH field 

First, the zonally averaged RH of January and July calculated 
from this scheme are compared with the climatology compiled by Thomp- 
son and Warren (1982), whose data are from London (1957) for 7~ 
22km, 0 - 75°N; Oort and Rasmussen (1971, hereafter denoted as OR) 
for 0 ~ 7km, 10°S 75°N and from Sasamori, et_ (1972) for 

0 ~ 22km, 10 ~ 90°S. 

Figures 4.5 and 4.6 show the comparison of the averaged RH of 
January and July respectively. It is clear that both cases have good 
agreement in low and mid-latitudes with a maximum at the tropics 
which coincides with the ITCZ, and minima at the subtropics which 
coincides with the descending arm of the Hadley cell. The seasonal 
variations of the pattern are also in phase with the position of the 
sun. Significant differences occur at the high latitudes, where 
the polar descending flow would result in lower _RH; the winter polar 
region is overestimated while the summer polar region is underesti- 
mated. A possible explanation of this difference for the Northern 
Hemisphere is that London (1957) only provided RH data up to the 
tropopause, so that the upper level data are unavailable in these 
regions and in the high latitude the tropopause is lower than 12km. 
Therefore in January, Thompson and Warren's (TW) calculation does 
not have relative humidity input between the tropopause and 12km in 
the high latitude region. As the result, their RH values at these 
regions are lower. However, in the July case the mixing ratio of the 
lower atmosphere in the Northern Hemisphere is much higher than that 
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of OR. Generally speaking the OR mixing ratio is higher than the 
the mixing ratio derived from the data used in this study for both 
months (see Tables 4.2 & 4.3). This contributes to higher relative 
humidity, especially when compared with dew point derived relative 
humidity as this study has done. Another source of error is that 
Harries equau m 4,2 does not account for seasonal variations. 

Fortunately, outgoing longwave radiation is less sensitive to 
averaged relative humidity when the surface temperature is lower, 
see Figure 2.1. In other words, the polar regions can tolerate 
more averted RH error than the tropics. Take the worst case, in 
July at 65°N, where the surface temperature is 12°C. If one varies 
relative humidity from 45% to 55%, clear sky outgoing longwave radi- 
ation only changes by about 5 W/n»2. Under cloudy sky conditions 
the error would be smaller. 

4.6.2 Zonal ly Averaged Clear Sky Outgoing Longwave Radiation 

In figures 4.7 and 4.8, the zonally averaged clear sky outgoing 
longwave (LW) radiation is compared with the calculation by Warren 
and Thompson (1983) for January and July. The excellent agreement 
between these calculations is obvious. The differences for most 
latitudes are smaller than 3W/m2. The zonally averaged surface tem- 
perature is also plotted on the same figures so that the relation 
between surface temperature and outgoing LW can be easily seen. 
Outgoing LW radiation versus RH is plotted in Figures 4.19 through 
4.21 for the July case. In mid and low latitudes where the surface 
temperature is high, it is very evident that the outgoing LW is 
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Figure 4.7 Zonal ly averaged clear sky OLWR and T s for January. 

The solid line is calculated from current study; dash 
line is from Tw and dash-dot line is T s . 
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Latitude 


Figure 4.8 same as figure 4.7 but for July. 
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anti correlated with jjtH. At high latitude, where the surface tem- 
perature is low, outgoing LW is insensitive to averaged RH but is 
nearly linear with the surface temperature as discussed in the pre- 
vious section, 

4,6,3 Global RM Distribution 

The world wide RH distributions are plotted in Figures 4.9 and 4.10 
for January and July respectively. Compared with temperature contours 
in Figures 4.11 and 4.12, the RH contours have much higher longitu- 
dinal variations. Still, for both months the subtropical subsidences, 
the so called dry tongue of the Hadley cell, can be easily identified 
with low RH centers, as well as the high RH centers along the tropics 
with the Congo, Indonesia and the Amazon being the most prominent 
three regions. In January over the ocean, there are high RH centers 
on the east sides of the Pacific and the Atlantic which result from 
the cold temperature of the California and Canaries currents. The 
low RH center at the northeastern Pacific is possibly due to the 
warm surface temperature brought by the northern bound Kuroshio cur- 
rent. Over the continents, the cold temperatures at Siberia and 
Greenland also contribute to high RH as the water vapor in the air 
is conserved. 

In July, the averaged R!H center along the subtropic subsidence 
belt of the Southern Hemisphere remains fairly well defined as in 
the January case. However, the low average RH center in Australia 
is more clear for the reason that in the Southern Hemisphere winter 
land-ocean temperature differences are more pronounced. The high 
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averaged RH centers along the intertropical convergence zone (ITCZ) 
all move northward. The region over the Amazon extends over the 
Caribbean sea. The low RH over the Sahara desert intensifies to cover 
the Middle East. The high RH center over Tibet and Central Asia 
weakens from a winter value of about 60% to 48% due to increasing 
air temperature and lack of influx of water vapor. The same applies 
to Siberia. 

4.6.4 Global Clear Sky Outgoing Longwave Radiation Distribution 

Figures 4.13 and 4.14 show the January and July clear sky outgoing 
longwave radiation field. It is easily seen that the winter hemi- 
sphere has a higher latitudinal LW radiation gradient than the summer 
hemisphere. The influence of the continents and of the ocean currents 
along the coasts can be easily identified. The highly longitudinal 
variations in the subtropics and tropics are remarkable. In this 
region, high LW centers are located in the central Pacific, Australia, 
and the southern Indian Ocean; and the low LW centers are located in 
Indonesia, the Congo and the Amazon, which coincides with the high 
averaged RH_ centers at those regions. These three centers move north- 
south according to the sun's position as well as change their inten- 
sity. Generally speaking, the high LW centers are located in the 
subtropics and low LW centers are located in the tropics. This agrees 
with the zonal ly averaged results discussed by Warren and Thompson 
(1983). 

An important feature of Figures 4.13 and 4.14, which cannot be 
seen in the zonally averaged plots is that along the tropics the 
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Figure 4.14 Clear sky OLWR map for July, in w/m 2 
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longitudinal variations are the strongest of the entire globe. Also, 
the high LW radiation at the central Pacific near the date line is 
as strong in the tropics as in the subtropics for both January and 
July and the intensity over the tropical Atlantic is at the same 
magnitude as in the subtropics as well. It is the three huge low LW 
centers of the Congo, the Amazon and Indonesia that contribute to 
the low zonal average outgoing LW at the tropics. It is somewhat 
misleading to say that the outgoing LW at the tropics is relatively 
low compared to that of the subtropics. 

Over the continents, reducing the air temperature will increase 
RH, so that low temperature with high relative humidity will further 
reduce low outgoing LW radiation. This phenomenon occurs in winter 
Siberia, Greenland and summer Central Asia, where the surface temper- 
ature of the high plateau is relatively low compared to adjacent 
regions, therefore a deep trough of LW contours occurs here. 

Minnis and Harrison (1984, hereafter denoted as MH) have studied 
a regional clear sky outgoing longwave radiation (OLWR) of November 
1978 using GOES data. The monthly mean of this region, from longitude 
30°W to 125°W, latitude 45°N to 45°S with 2.5°x2.5° resolution is 
shown in Figure 4.15. For comparison, our calculation of climatolo- 
gical November clear sky outgoing LW flux is shown in Figure 4.16. 
Regardless of the coarser resolution of our coordinates and the dif- 
ferent averaging time, there is good agreement between the two results. 
The locations of the low OLWR flux center at northwest South America, 
the high OLWR flux centers off-shore of Brazil, the southeast Pacific 
and Caribbean sea, and a shallow trough over continental United States 
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are all well simulated. Also the Intensity of OLWR flux over the 
oceans are In good agreement. However, over Brazil our calculation 
is about lOW/m 2 overestimated. This might suggest that averaged £H 
Is underestimated at this region, or the surface temperature of 
November 1978 is below normal at this region. 

As one compares the OLWR field, Figure 4.16, with temperature and 
averaged RH fields, Figures 4.17 and 4.18, it is clear that averaged 
RH highly changes the OLWR flux pattern. This result indicates that 
water vapor is one of the dominant factors of outgoing LW flux. 
Even without the presence of cloud, OLWR field is modulated by 
water vapor distributions. 

The simple Budykc-Sel lers type parameterization of OLWR calcula- 
tion is often used: 

F* A+BT 4.12 

where F is OLWR flux, T is surface temperature, and A and B are 
constants. The above result implies that this equation significantly 
biases the radiation field by omitting water vapor. 

4.7 Sensitivities of the Interpolation Coefficients 

As pointed out in section 4.3, some coefficients in equation 4.2 
and 4.4 are rather arbitrarily chosen. These might have substantial 
impact on the outgoing LW radiation. In this section, we present a 
sensitivity analysis of coefficients to show how the outgoing LW 
responds and verify that the outgoing LW radiation has minimal sen- 
sitivities to the changes of these coefficients. For convenience, 
these two equations are repeated here, and only the results of the 
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July cases are plotted. The January cases have similar results. 

OlOOnto ( J 5 -3xl°- 6 /(1+ ( J -20 ) 2 /i2 8 ) 4.2 

n(J)=0.86x(l+(J-19) 2 /1500) 4.4 

Case A 

The coefficient 0.86 in equation 4.5 is perturbed to 0.80 and 
0.90. The results of the July case are plotted in Figure 4.19. As 
the figure shows, changing the coefficient 0.86 has more influence 
on the tropics and less toward the poles. Changing the coefficient 
from 0.86 to 0.80 increases the RH at the equator by 9% and reduces 
the clear sky outgoing LW by 6 W/m 2 . Changing it from 0.86 to 0.90 
reduces the RH at the tropics by 5% and increases the clear sky 
outgoing LW by 2.5W/m 2 . These differences are smaller than those 
from di fferent radiative transfer models, where 8W/m 2 variation 
on OLWR are not un'-ommon (Luther, 1983). The original number used 
by Briegleb and Ramanathan (1982) is 0.85. The differences between 
their work and this study art negligible. 

Case B 

Perturbing the coefficient 1500 in equation 4.5 from 1000 to 
2000 only slightly influences the RH and outgoing LW at mid latitudes, 
as shown in Figure 4.20. This is the least influential coefficient. 

Case C 


In Harries' water vapor mixing ratio formula, equation 4.2, the 
coefficient 3x10“^ stands for the average mixing ratio measured by 


S 80 60 40 20 0 20 40 60 80 N 
Latitude 

Figure 4J9 Sensitivity of _RH and clear sky OlWR to the interpolation 
parameter, Case~A. XIV is the undisturbed case; II changes 
0.86 in equ. 4.5 to 0.80 and III changes 0.86 to 0.90. 
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Figure 4.20 Same as figure 4.19, but for Case B. XIV is the undisturbed 
case; VI changes 1500 in equ. 4.15 to 2000 and VII rS anges 
1500 to 1000. 
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various experiments at lOOmb, see Figure 4.1. This coefficient has 
a linear response or the averaged Rji for all the latitudes. For 
example, at the tropics, changing the coefficient 3x10“® to 6x10“®, 
see Figure 4.21, would increase the averaged J1H from 54% to 60% and 
reduce clear sky outgoing longwave radiation by approximate 4W/rr»2. 
Changing the coefficient 3x10“® to 1x10“® would reduce the averaged 
RH from 54% to 47% and increases clear sky outgoing longwave radi- 
ation by about 5W/m2. Therefore, a six-fold change of mixing 
ratio, from 1x10"® to 6x10“®, would only change clear sky outgoing 
longwave flux by 9 W/m^, only about a 4% variation. Thus the 
longwave flux calculation is insensitive to this perturbation. 

In conclusion, we adapted all the coefficients chosen by other 
investigators except to slightly tune the coefficient C.85 in equation 
4.5 to 0.86. We found that the averaged _RH profile is in good agree- 
ment with Thompson and Warren's in the low and mid latitudes, and 
also found ver^ good agreement in zonal averaged clear sky 0LWR 
profile for all latitudes. The changes in these coefficients has 
little effect on clear sky 0LWR which would imply even less error 
for cloudy sky 0LWR. '' 1 oudy sky 0LWR will be discussed in Chapter 6. 
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Figure 4.21 Same as figure 4.13, but for Case C. XIV is the undisturbed 
case; XI changes 3x10"^ in equ. 4.2 to 6x10“® and X changes 
3xl0“ 6 to lxlO" 6 . 


CHAPTER 5 


CLOUD AND CLOUD TOP TEMPERATURE 

This Chapter is devoted to two subjects: cloud cover end cloud 
top temperature T c . For the cloud cover, from section 5,1 through 
5.4, the impact of cloud on climate and cloud observation is first 
reviewed. Then, the Sherr et^ aj_. (1968) cloud cover model, which is 
incorporated into the outgoing longwave radiation (OLWR) calcula- 
tion, is discussed. The distribution of cloud cover follows. For the 
cloud top temperature T c , in section 5.5 the techniques of the mea- 
surement are first reviewed, then a set of T c data is compiled and 
compared with a region with satellite observation. A brief summary 
follows in section 5.6. 

5.1 Cloud and Climate 

Earth OLWR, unlike idealistic clear sky OLWR as discussed in the 
previous chapter, is highly modulated by the dynamics nd thermo- 
dynamic processes of clouds. The presence of most clouds traps the 
higher longwave (LW) radiation from the surface and emits its own 
lower LW radiation to space. This process conserves earth OLWR, and 
is often referred to as the cloud greenhouse effect. Wit.i a theore- 
tical model, Kuhn (1978) showed that outgoing longwave flux (OLWF) 
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depends on cloud height, e.g. the greenhouse effect intensifies as 
the elevation of the cloud deck increase. However, in polar regions 
where the temperature inversion is strong and water vapor density is 
low, inserting a cloud deck would increase OLWF substantially. 
More recent observations by Minnis and Harrison (1983) also showed 
high anti -correlation between diurnal cloud cover and OLWF for a 
few locations in the mid latitude of the Southern Hemisphere. The 
multilayered clouds complicate this problem, especially high thin 
cirrus type clouds with less opacity. For example, it was pointed 
out by Hetman (1981) th*t sophisticated GCM's underestimated OLWR 
because higher LW radiation from lower warmer cloud can penetrate 
througn overlying cirrus. 

In addition to its effect on the LW radiation, cloud also has 
very strong influence on solar short wave (SW) radiation. Its pre- 
sence reduces insolation on the surface and its albedo, which is 
higher than the surface, increases tne reflection of SW to space, 
which is often referred to as cloud albedo effect. 

Global surface albedo is only about 13%, while cloud albedo can 
range from 20% to as high as 90%. This effect competes with the 
greenhouse effect, so that the degree of the cloud's net ; mpact on 
the climate, cloud feedback effect, is still controversial . Schneider 
(1972) first used an atmospheric model to analyze the relative mag- 
nitude of each compo n ent. He defined a parameter 
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where A is absorbed solar flux, F is OIWF, and C is cloud cover. 
When 5 is negative, the albedo effect dominates, and when 6 is 
positive, the greenhouse effect dominates. He found that, with the 
constraints of fixing other atmospheric factors, the effect of in- 
creasing global average cloud cover is to decrease the global average 
surface temperature due to higher albedo effect. However, Cess and 
Ramanathan (1978) stressed that any discussion of the radiative 
effect of changing global cloud amount or atmospheric transparency 
must consider the possibility of change in the distribution and the 
type of cloud as well, which is not well understood yet. One of the 
most difficult problems involved is the the lack of an unique defi- 
nition of cloud cover. Because cl -ud has different characterise cs 
for different wavelengths, effective cloud cover for incoming solar 
radiation might differ significantly from that of outgoing LW radia- 
tion. Ohring and Clapp (1980) and Hartman and Short (1980) used 
OLWF versus albedo instead of using cloud cover directly to determine 

the net energy budget. 
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where a is albedo, Q is one fourth of the solar constant. This 
method, however, did not eliminate other atmospheric variations, 
such as water vapor content, lapse rate, etc.. Cess, et aj_. (1982) 
also d : cussed that using different satellite data would change the 
results substantially due to different data sampling technique. 
Although the effect of the cloud feedback on net radiation balance 
is still uncertain (Hartman, 1983), there seems to be a good agreement 
with the result of Schneider (1972) that cIouj albedo effect domi- 
nates the greenhouse effect except for very high cirrus type cloud. 

Warren, et al_. (1983) used 35 years of ship observations to 
study cloud cover over ocean area, and found that there are statis- 
tically significant trends for different types of cloud over different 
areas, e.g. 11 per 10 year decrease in cumulus over most ocean 
areas of both hemispheres and a similar increase in stratus and 
stratocumulus over all ocean areas except the northwest Atlantic 
during all seasons and the eastern South Tropical Atlantic during 
Southern Hemispheric summer. It is intriguing to ask how this 
trend changes the climate or why this cloudiness trend appears in 
the data. But, the first important problem involved here is cloud 
observation. 

5 .2 Cloudiness Observation 

Tt is redundunt to stress how subjective cloud observations, both 
cloud classification and cloud cover, are. A simple example is that 
cloud cover reported in early morning and late afternoon tend to be 
overestimated because high solar zenith angle will enhance the perce- 
ption of the lateral area of cloud. Another problem on surface 
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observation is the lack of stations over many remote regions, although 
ship observations help to fill in very sparse data over ocean. To 
have a global surface cloud observation network is very impractical. 
Despite such difficulties, London (1957) compiled a seasonal clima- 
tology of total cloud cover distribution for the Northern Hemisphere. 
It is still one of the most referred works, regardless of their 
accuracy. Later a similar work for the Southern Hemisphere was done 
by van Loon (1972). 

Satellites provide an easier way for sampling data. Early polar 
orbiting TIROS satellites had limited ground coverage and limited 
sampling times. Clapp (1964) utilized this satellite photography 
and performed some nephanalyses of seasonal global cloud cover from 
March 1962 to February 1963. This study showed that global cloud 
has some characteristic patterns over seasons and concluded that 
the cloud distribution on large-scale must posses a certain stability 
or repetitiveness over the entire seasons, even in the regions of 
migratory cyclones and anticyclones. Clapp did not quantitize the 
result, such that it can not be used to verify London's (1957) 
previous work; nevertheless it did verify the feasibility of classi- 
fication of cloud cover from space and its agreement with surface 
observation for zonal ly averaged results. Another work by Kornfield 
and Hasler (1969) used a multiple exposure technique for developing 
monthly averaged cloud cover pictures and then performing nephanaly- 
ses. Ths study clearly identified the ITCZ and some quasi -stationary 
waves. Still both Clapp's and Kornfield and Hasler's work were 
very dependent on subjective judgement of an experienced analyst. 
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A number of more objective methods have been proposed. Instead of 
using photography, different channels of radiation band are used tc 
reconstruct cloud pictures, such as was done by Shenk and Solomonsor 
(1972). Their threshold method for using visible channel brightness 
from a high resolution scanning radiometer can produce good results, 
if the field of view is completely filled with cloud or completely 
cloud free. However, if the assumed condition does not prevail, the 
potential error could be as high as 75%. Another method, the spatial 
coherence method (Coakley and Bretherton, 1982) using advanced very 
high resolution radiometer (AVHRR) information, can determine single- 
layered cloud cover or multilayered cloud cover extended over moderate 
layer regions. Chahine et al_. (1983) used multichannel high resolu- 
tion infrared scanner (HIRS) and microwave sounding unit (MSU). This 
method can determine the effective cloud cover under the assumption 
that the cloud is single-layered and non-reflecti ve. 

An effort, the International Satellite Cloud Climatology Project 
(ISCCP, Schiffer and Rossow, 1983) has been initiated for organizing 
che data collection and processing to determine temporal and spatial 
variability of the amount, type and radiative properites of clouds. 
This study will not only provide a global data set for clouds, but 
also lead toward a better understanding and enhan^r the ability of 
cl imate modeling. 

5.3 Sherr et al. Cloud Model 
5.3.1 The Model 

Sherr et al. (1968) compiled an empirical global cloud cover mo- 













del, originally intended for simulation of the effects of cloud 
cover on proposed earth-viewing space missions. The surface of 
earth is divided into 29 cloud climatological regions as shown ;* 
Fig. 5.1, Within each region, cloud is assumed homogeneous. This 

feature provides the interface with the outgoing LW parameterization 
of Thompson and Warren (1982) used in this investigation, which 

requires a single deck effective cloud cover as a predictor (see 

Chapter 2). We adapted the model with modifications suitable for 

this study. A general description of the nature of the cloud in 

each cloud climatic region is given in Table 5.1. 

In constructing the model, Sherr et (1968), the data from 

satellites were used to identify the regions where the cloudiness 
is homogeneous. Ground-based observations were used to construct 

cloud cover frequencies. Data from as man w areas as possible were 
assembled and from these a representative station was identified. 
Data from this station were usually based on at least 10 years of 
observations. 


5.3.2 Other Studies with Sherr et al . Model 

Based on the results of Sherr et _al_. , Falls (1974) statistically 
analyzed the cloud cover frequency and concluded that world cloud 
cover behaves as a Beta probability distribution, which can describe 
the cloud cover field well with relevent values of two parameters. 

Based on the surface classification of Sherr et aj_. model, 
Henderson-Sellers (1978) studied the surface type and its effect on 
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> 70 

Synoptic 

Scale 

Small 


< 
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Or 


Table 5.1 continued. 


1 

2 

3 

4 

S 

6 

7 

8 

9 

14 

High Latitude 

Northern 

Moderate 

> 80 

■W 70 

Synoptic 

Small 

.... 

- Ocean 

Hoal sphere 




Scale 



IS 

Mar 

Northern 

Small 

n. 80 

n. 80 

Synoptic 

Small 

.... 


Heel sphere 




Scale 



16 

Tropical - 

North of 

Extreme 

> 70 

< 40 

Convtctlve 

Large 

1600 


Seasonal 

Change 

Region 03 







17 

Tropical - 

Northern 

Moderate 

^ 50 

< 30 

Cbnvectlve 

Large 

1600 


Clear 

Hemi sphere 








Winter 

Near Region 16 







18 

Ned iter- 

Northern 

Extreme 

n. 30 

.. 

Convective 

Small 



ranean 

Hemisphere 

Europe, 



' 80 

Synoptic 

Small 

• ••• 



Western North 
America 




Scale 



19 

Sub Tropical 

Northern 

Moderate 

•v. 50 

-- 

Convective 

Large 

1600 


Hemisphere 


— 

> 60 

Synoptic 

Small 





n- 30N 




Scale 



20 

Sub Tropical 

Northern 

Moderate 

n. 50 

-- 

Convective 

Small 

— 


- Ocean 

Hemisphere 


— 

» 60 

Synoptic 

Small 

— • 



■v 3 ON 




Scale 



21 

Tropical - 

South of 

Moderate 

-u so 

> 60 

Convective 

Large 

1600 


Cloudy 

Scanner 

Region 03 







22 

Mid Latitude 

Southern 

Moderate 

■x. 70 

■*. 60 

Synoptic 

Small 

.... 


Ocean 

Hemlspnero 




Scale 



23 

High Latitude 

Scuinern 

Moderate 

•>. 70 

> 80 

Synoptic 

Small 

— 


• Ocean 

Hemisphere 




Scale 



24 

Polar 

Southern 

Small 

•v. 60 

i 60 

Synoptic 

Small 

— 



Hemisphere 




Scale 



2S 

Tropical - 

South of 

Extreme 

< 40 

» 70 

Convective 

Large 

1600 


Seasonal 

Change 

Region 03 







26 

Tropical - 

South of 

Moderate 

« 30 

% 50 

Convective 

Large 

1600 


Clear 

Region 25: 








Winter 

Africa, 

Australia 







27 

Hedlter. 

Southern 

Extreme 


30 

Convective 

Small 

— 


ranean 

Hemisphere 
Austral ia , 


n 60 


Synoptic 

Small 

.... 



Chile 




Scale 



26 

Sub Tropical 

Southern 

Moderate 

.. 

< 50 

Convective 

Large 

1600 


Land 

Hemisphere 


•V. 50 

-- 

Synoptic 

Small 




•v 30S 




Scale 



29 

Sub Tropical 

Southern 

Moderate 

.. 

«W 50 

Convective 

Small 

— 


• Ocean 

Hemisphere 


> 60 

— 

Synoptic 

Small 

---- 



% 3 OS 




Scale 
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cloud cover. She has also found the Beta distribution useful as a 
model for the cloud cover frequency. Her results also imply that 
cloud frequency curves over land in high latitudes, :.n d-1 at i tudes 
and subtropical latitudes are similar, while the frequency distribu- 
tions of cloud cover differ (1) over land and ocean at the same lati- 
tude and (2) over ocean areas in different latitude zones. 

Combining with 45 months of satellite radiation measurements, Bean 
and Somerville (1981) derived a Beta distribution world-wide cloud 
cover model . This model can describe cloud fields very concisely. 
However, the parameters were derived from fairly small size sample, 
such that statistical confidence is not achieved. Also, utilizing 
the model to other time period other than that the parameters were 
derived might cause significant error. 

5.3.3 Modification of Sherr et al. Mode. 

A version of the Sherr et aj_. cloud model used by Bartman (1980} 
in his albedo model is incorporated into this study. Since the cloud 
climatic region areas of Sherr et_ al_. did not begin and end on the 
10°xl0° latitude-longitude lines of our models, Bartman stretched cr 
contracted slightly the map areas in order to fit the model grid 
system. Table 5.2 shows the cloud climate regions used for each 
10°xl0° area in Bartman 's model and this study. 

As the result of this process, one would notice that along the 
latitude 65°S on table 5.2, the cloud typ . is alternating between 
type 23 and type 24. On a contour plot of cloudiness, this would 
produce unrealistic wavy curves along latitudes 55° to 75°. Thus, 
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instead of directly using the cloud cover at this region generated 
by the classification, we used equations 5.3 and 5.4 to smoothe the 
distribution portrayed in fig. 5.1. 

f ( J»16)*f (j , 17 )*2/3+f ( j ,15)*l/3 1 < j < 7 or 19 < j < 36 5.3 

f (j *16)*f (j ,17 )*l/3+f ( j ,15)*2/3 8 < j < 18 5.4 

where j is a longitude index such that longitude is J* 10°-5°, and K 
is a latitude index such that latitude is 85°-K*10°. 

Monthly cloud cover used in this study is calculated as follows. 
Sherr et_ al_. world-wide cloud cover distributions are given in five 
cloud cover categories, listed in table 5.3, for each of the 29 cloud 
climatological regions. Frequency distribution of the five possible 
cloud amount categories are given for eight local times at 0100, 
0400, 0700,..., 2200. From the frequency distribution data and 
cloud cover categories, the percent cloud cover was calculated for 
each cloud climate region as a function of the eight different local 
times of the typical day for each month of the year. The monthly 
cloud cover is obtained by averaging the cloud covers of the eight 
local times. 

A calculation of the global annual percent cloud cover from the 
above data gave a value of .60, which is much larger than the values 
.50 to .55 normally used, e.g. the global cloud cover of Hoyt (1976) 
was .532; C~ss and Ramanathan (1978) suggested annual global cloud 
cover is between .45 to .55. Ramanathan (1976) used .45 as the effec- 
tive total cloudiness for his radiative convective model. 
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Category 


Cloud Cover in Tenths 


1 

2 
3 
k 
5 


0 

1,2,3 

4,5 

6,7, 8,9 
10 


Table 5.3 Cloud Cover Category 
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A multiplying factor 6LCLC/60 has thus been incorporated into 
the cloud cover calculation. This factor is not only used to adjust 
the difference between the Sherr et al_. model and the nominal value, 
but also accounts for the emissivity of the cloud, such that GLCLC/60x 
A c equals effective cloudiness, where A c is cloud cover from the 
Sherr et £1_. model. In the current study, GLCLC is set to 45 as by 
Bartman (1980, 1981). 

5.4 Cloud Cover Distribution 
5.4.1 Zonally Averaged Cloud Cover Distribution 

Figures 5.2 and 5.3 show the zonally averaged cloud cover from 
(1) the modified Sherr' s cloud model without the factor GLCLC and 
from (?) London and van Loon for January and July respectively. 

Both months have general agreement between (1) and (2), i.e. maxi- 
mum cloud cover in the tropics and middle latitudes and minimum 
cloud cover in the subtropics and polar regions. The high correlation 
between cloud cover and uppper air relative humidity can also be 
well identified (see Figures 5.4 and 5.5, taken from Telegadas and 
London, 1954 ) . 

However, the differences of the magnitudes between (1) and (2) 
are also very significant. Poleward greater difference could be 
possibly due to less observation frequency or higher human perception 
bias from larger zenith angle, or both. The verdict of better 
cloud cover climatology will not be issued until either Warren et al. 
or ISCCP results are revealed. 


-90 -60 -30 0 30 60 90 

Latitude 

f' a ure 5.2 Comparison of the zonally averaged cloud cover between 
the climatology and Sherr £t £l_. (1968), January. 
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Cloud cover 




10 


to 


to 


to 


to ac *o 


0 

Figure 5. 


LATITUDE *N — 

Averaged latitudinal distribution of the total cloud 
cover and relative humidity at the 850, 700 and 
500mb-- -wi nter after Telegadas and London (1954). 
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5.4.2 Global Distribution of Cloud Cover 

Figure 5.6 and 5.7 are cloud cover contours from the data of Sherr 
et al_. (1968) for January and July respectively. 

In the GlwO. l :loud distribution, maximum cloud cover often 
occurs at the *TCZ, and the mid-latitude westerly zones where there 
is a warm moist surface region or large baroclinicity, etc., while 
minimum cloud cover occurs at the subtropical subsidence, high pre- 
ssure centers, and polar subsidence regions. Other mechanisms, such 
as monsoonal circulation, also generate extreme cloud cover. 

Tables 5.4 and 5.5 summarize the causes and the regions of extreme 
cloud cover for January and July respectively. Although for some re- 
gion there might be more than one reason for the causes of extreme 
cloud cover, e.g. ITCZ coincides with moist surface, only the dominant 
reason is mentioned. 

5.5 Cloud Top Temperature 
5.5,1 Definition and Measurement 

As pointed out both in Chapter 2 and in the previous section, the 
cloud cover used in this study is defined as the effective cloud cover, 
with emissivity equal to one. Thus, the definition of cloud top 
temperature refers to the temperature at the top of the effective 
cloud. Since not all of the cloud is black, the virtual cloud top 
might be lower than its physical location in order that emission at 
warmer level would compensate the smaller emissivity. However, in 
practice the exact location of the effective cloud top can not be 
defined using conventional radiosonde, neither can the cloud top 
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January 

Causes 

Locations 

maximum 

cloudiness 

ITCZ 

equatorial low 

1. Maximum cloudiness along 5*S In the 
Indian Ocean continues Into Western 
Pacific. 

2. Eastern Pacific & Atlantic In the tropics 
which coincide with the strongest trade 
winds and cold sea surface temperature 
except Indian Ocean. 

Baroclir.icity 

1. In the Northern Hemisphere, the prevailing 
westerlies and their cyclonic storms are 
located at North America, Nortwest Europe. 
Northwest coast of North America, Alaska 
and Far East. 

2. In the Southern Hemisphere are at 
circumpolar zone at subpolar latitudes. 

warm moist 
surface 

1. The amazon and Congo Basin where tropical 
forest located. 

2. Indonesia & tropical coast of Australia 
where warm sea surface temperatures 
located. This Is also the region where 
the ascending arm of Walker circulation 
is located. 

minimum 

cloudiness 

subtropical 
high pressure 

Eastern and southwestern section of 
subtropical high pressure in the Northern 
Hemisphere, eastern and northeastern 
section In the Southern Hemispheric high, 
Sahara, West Africa desert Southwest 
United States and Southwest South America. 
In the Southern Hemisphere the high 
pressure spread over almost all longitudes 
except those of the Andes mountains. 

monsoonal 

circulation 

India under northeasterly winds. 

Polar region 

North and South poles. However, it is 
difficult to confirm the observation 
during polar night. In the Northeast 
Siberia, Arctic and northern America, the 
cloudiness is highly dependent on the 
amount of open water. 


Table 5.4 Summary of extreme cloud for January. 





















"TuTy 

Causes 

Locations 

1 


1TC2 

equatorial low 


Near equator in the Indian Ocean and an 
extentior. of cloudy band southeastward 
from New Guinea. 

maximum 

cloudiness 


2. 

Northern and Central jouth America, 
Central Africa and equatorial Pacific. 

Baroclinicity 


North displacement of Bermuda and Pacific 
high pressure centers: Principal cyclonic 

activities over Europe and North America 
is confined to the northern regions off 
Alaska, northern Europe, central and 
northern United States. 


monsoonal 

circulation 

India 

minimum 

cloudiness 

subtropical 
high pressure 

1. 

Northern Africa, southwest Asia, central 
Pacific and southwest United States where 
are on the eastern slopes of the high 
pressure waters. 



2. 

South Pacific and West shore of South 
America, Desert Australia. 



Subpolar circumcircle over the Southern 
Hemisphere has less cloudiness than 
January. 


Table 5.5 Summary of extreme cloud cover for July 
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temperature be directly measured. 

Cloud top temperature can be indirectly measured from satellite. 
As in the common single channel threshold method, cloud top temper- 
ature is the infrared equivalent blackbody temperature of cloud, 
which is measured by using a window band radiometer. Minnis and 
Harrison (1984) have analyzed a set of GOES data using this techinque. 

A recent proposal by Chahine et aj_. (1983) presented a method 
using multichannel HIRS to retrieve both effective cloud cover and 
cloud top temperature independently. Under the assumption that the 
reflectivity of cloud in the i'R region is zero and that the emissivity 
of the cloud at different infrared wavelengths are the same, the 
pressure level for the black cloud top location can be retrieved. 
The cloud top temperature is then assumed equal to its ambient tem- 
perature. The corresponding temperature sounding is retrieved using 
other channels. Although this method has been tested on NOAA weather 
satellite data, a global monthly average climatology is not yet 
avai lable. 

Since cloud top temperature is the most indefinite parameter in 
calculating outgoing LW radiation, it is taken as an adjusting factor, 
as suggested by Thompson and Warren (1982), to match the current 
climate. 

5.5.2 Compilation of Cloud Top Height and Cloud Top Temperature 

According to the cloud region classification of Sherr et _al_. 
(1968), section 5.3, the cloud cover of each of the 29 cloud climato- 
logical regions is homogeneous. Table 5.1 also gives the general 
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description of the cloud type for each region. Since each cloud 
type has its prevailing cloud top height, one can assume that in 
each region all the cloud top Is at the same altitude, and then take 
the cloud top ambient temperature as the cloud top temperature T c 
as was done by Hummel and Kuhn (1983) and Chahine et aj_. (1983). 

To avoid generating more uncertainty, cloud top height is only 
assigned to the mandatory pressure levels, where the temperature 
field Is written on the NCAR tape (see Chapter 3). Since there is 
scarcely any cloud sustained at lOOmb, no cloud is assigned to this 
level. Also, very low level clouds are assigned to the midpoint 

between the surface and 850mb. Thus each cloud top is at one of the 

following levels: near surface, 850, 700, 500, 300 and 200mb. Sea- 
sonal variation of cloud top height is cc'^idered as well. Table 

5.7 shows the cloud top level of each clou 1 type for each month. 
Determination of the cloud top level for each cloud region will be 
discussed in the following paragraphs. The temperature at this 
pressure level is taken as T c , except for the near surface cloud 
where the mean temperature between the surface and 850mb is used. 

An atlas of monthly 0LWR is compiled as the reference for determ- 
ining the proper cloud top height for each cloud type given in Table 
5.1. Data used are from the scanning radiometer (SR) on board of 
N0AA satellites measured from June 1974 to February 1978 (Abel and 
Gruber, 1979). These data were further processed into 10°xl0° regions 
by Vonder Haar et aj_. (1982). This is the longest time series 
of satellite measurement from one single instrument. Data of the 
same month from different years are averaged to generate the 'monthly 
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average' . 

To determine the cloud top level, a rough guess was first made 
by using the description of Sherr et j*l_. (1968) for each cloud type. 
Table 5.1, and choosing the possible corresponding cloud height level 
from Table 5.6 from the International Cloud Atlas (1975). Thus, 
T c can be found, then the cloudy sky OtWF is calculated. This 
result is compared with the reference atlas, cloud region by cloud 
region, and the cloud top level can be adjusted to the upper level 
or the lower level In order to minimize the difference between the 
reference atlas and calculated OLWF. This procedure is iterated 
until the difference is minimal for each cloud region, and the best 
cloud top height is thus defined. Table 5.7 tabulates the pressure 
level where the cloud top is located, and is the result of this 
procedure. 

Examining Table 5.7, one finds that cloud climatological regions 
have seasonal variation of cloud top height. The most drastic one is 
Region 25, located at the southern tropics over South America, Africa 
and Indian Ocean, where the cloud top level changes from 200mb in 
January to 850mb In July. The range of the seasonal variation of 
cloud top height is about 200mb to 300mb for most other cloud regions. 

There are as well many regions where their cloud top height 
stays fairly constant throughout the year, such as region 3 at the 
tropical Pacific and the northern part of the South America; Regions 19 
and 29, where subtropical oceans are located; Regions 14 and 23 of 
the high latitude ocean; Region 11, mid-latitude land in the Northern 
Hemisphere; Region 23, mid-latitude ocean in the Southern Hemisphere 
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and Region 28, the subtropical land at the southeast parts of both Aus- 
tralia and South America. 

Figures 5.8 and 5.9 show the contour of cloud top temperature for 
January and July respecti vely. It is clear that the cloud region 
classification of Sherr et_ aU (1968) highly modulates this result. 
In January, the main low cloud top temperature regions are located 
at the west coast of the United States, Siberia, equatorial Pacific 
and Atlantic, northern South America and African Congo. High cloud 
top temperature regions are Australia, west coast of Central America, 
west coast of South Africa, monsoonal India and West Europe. In 
July the main low cloud top temperature regions are located at the 
eastern United States, Caribbean Sea, southeast Pacific, East China, 
and monsoonal India. High cloud temperature regions are located at 
North and South Africa, southern Indian Ocean, west coast of the 
United States and eastern South America. 

As mentioned in the previous section on cloud top level, the 
regions with strong seasonal variation are South America south of the 
equator, monsoonal India, and South Africa. Other regions also have 
seasonal variation, but nevertheless, are more moderate. 

5.5.3 Comparison with Sate! 1 ite Observation 

As was done for the clear sky OLWR case in section 4.6, the com- 
piled T c is compared with satellite data analyzed by Minnis and 
Harrison (1984, hereafter denoted as MH) for the region from Longi- 
tude, 30°W to 125°W, latitude 45°N to 45°S for November. These results 
are shown in Figure 5.10 and 5.11. 
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Figure 5.10 GOES regional T c for November. 1978 after Minnis and 
Harrison (1984) , in °K. 








The locations of the region of the low T c center at central South 
America and the adjacent high T c regions both to the east over Atlan- 
tic and to the west over Pacific are well simulated. However, the 
magnitudes differ. In our compilation the maximum temperature Is 
294°K over the Atlantic versus 270°K of MH; the minimum Is 240°K 
for the Central South America, while MH has 260°K. At the Atlantic 

coast of South America, where the compiled data are higher, the 

reason for the difference might be that the NOAA data has higher 

OLWR at this region, 282 W/m^ vs. 270W/m2 for MH. Tne other possibi- 
lity Is that MH measured T c using a window channel which could be 
contaminated by the tropical water vapor absorption, thus underesti- 
mated. Over Central America, the compilation data is lower. This 
could possibly be explained as follows: our clear sky OLWR case in 
this region is about lGW/m^ higher than that of MH, which means 

the height weighted mean relative humidity _RH could have been under- 
estimated at this region. Therefore, to fit the normal cloudy sky 
OLWR from NOAA data, cloud must be adjusted to higher level in order 
to offset the low RjH, thus resulting in lower T c . Still, this 
magnitude is in the reasonable range. Comparing with other observa- 
tions, Curran and Wu (1982) used Sky lab-observed data and found T c 
can be as low as 226°K in New Mexico area. This somewhat justifies 
our T c compilation. 

5.6 Summary 

Clouds have two competing effects on climate: cloud albedo effect 
and greenhouse effect. Most observation shows that cloud albedo 
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effect dominates. However* the result Is sensitive to the analysis 
techniques and the data used. Since these cloud effects are highly 
dependent on the cloud type and cloud height, this problem can be 
very localized. Apparently, further study on this subject Is needed 
for Its great Importance on the earth radiation budget. On the 
large scale, satellite observations show that clouds possess certain 
stability or repetitiveness over an entire season, which render the 
feasibility of constructing cloud models, such as Sherr jet a]_. 
(1968) Incorporated in this study. Some techniques of cloud obser- 
vations are reviewed. 

The Sherr et a]_. (1968) cloud model classifies the globe Into 29 
homogeneous cloud cover regions. A set of monthly cloud cover maps is 
derived from this model with an adjusting parameter GLCLC introduced 
by Bartman (1980) for fitting the nominal global annual averaged 
cloud cover. A value of 45 Is selected for current study. The zonal 
averaged results are compared with London (1957) and van Loon (1972) 
with fair agreement. It Is expected that the future result of ISCCP 
and Warren et_ al_. would greatly Improve our understanding of cloudi- 
ness. 

A set of cloud top temperature maps are compiled using the cloud 
region classification of Sherr et aK (1968) and with the aid of the 
LW radiation measurement by NOAA scanning radiometer. A November 
regional result is compared with Minnis and Harrison (1984) with good 
agreement. 


CHAPTER 6 


CLOUDY SKY OUTGOING LONGWAVE RADIATION 

This Chapter presents the results of the cloudy sky outgoing 
longwave radiation (OLWR) calculated from the scheme developed in this 
study as described in the previous chapters. In section 6.1, the glo- 
bal averaged OLWR variation with season is discussed. In section 
6.2, the zonally averaged OLWR for January and July, as well as the 
temporal zonal variation are discussed. In section 6.3, the global 
distribution of OLWR for January and July are discussed in detail 
from region to region. 

6.1 Global Average OLWR 

Global monthly average cloudy sky OLWR calculated from the model 
(CD), clear sky OLWR (CR) and reference OLWR compiled from NOAA 
scanner radiometer (NOA, see section 5.5.2), are tabulated in Table 
6.1. A summary of other parameters are also included. 

The global average manifests the general performance of the model 
and also shows the mean state of the earth climate. There is good 
agreement between CD and NOA, with CD very slightly overestimated 
with an annual average of 245.6 and 244.7 W/m? for CD and NOA respec- 
tively. The differences are smaller than 1.5 W/m2, about 0.6%, 
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except for November, which differs by about 1%. The difference of 
the annual averages is about 0.9w/m2. As both cloudiness and mean 
relative humidity hold very steady throughout the year, at about 45% 
(after multiplied by GLCLC/60) for CD and and 47% for RH, the major 
parameters contributing to the seasonal variation are the change of 
surface temperature (T s ) and the temperature difference between 
the cloud top and the surface (T s -T c ). From Table 6.1 it is clear 
that both of t iese factors contribute to the maximum OLWR at the 
northern hemispheric summer, June through August. However, these 
two parameters are also slightly out of phase; T s has a maximum at 
July, while T s -T c has a minimum in May, with the result that June 
has the highest OLWR. Comparing CR with CD, on the average, clouds 
contribute about 10w/m2, 0 r 4%, in reducing OLWR. 

To illustrate the seasonal variation, CP and NOA are plotted on 
the climatology by Campbell and Vonder Haar (1980, denoted as CV) 
on Figure 6.1. The difference between CV and NOA will be discussed 
later. The seasonal cycle of OLWR is in phase with the season of the 
Northern Hemisphere with an amplitude of about 5 W/m^. The two prime 
factors affecting this seasonal variation are (1) Earth-Sun distance 
and (2) different amounts of land-sea distribution in the Northern 
Hemisphere and the Southern Hemisphere. Since the Earth-Sun distance 
is shortest in December, in the northern hemispheric winter, the 
factor of the land-sea distribution must dominate in order for the 
OLWR variation seasonal to be in phase with the northern hemispheric 
seasons. 

Although seasonal OLWR has a maximum in the Northern Hemisphere 
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Figure 6.1 Experiment comparison of monthly global average OLWR. 

Experiment 1 through 7 and 0 are from Campbell and 
Vonder Haar (1980). 


Ill 


summer and minimum in the Northern Hemispheric winter* the exact 
timing of the occurence for these extremes varies very substantially 
from year to year. A recent observation of Nimbus 6 and Nimbus 7 
from July 1975 to October 1980 by Bess (1984, personal communication) 
is shown on Figure 6.2. Disregarding an instrument calibration shift 
and missing data, the minima can occur from November to February and 
the maxima from July to August. NOA and CV, however, have maxima at 
June. The significance of this minor variation of OLWR to the climate 
can be interpreted in terms of net radiation, which is the difference 
between the two large values of 01WR and the absorbed radiation. An 
1 w/m2 net radiation flux would change the temperature of a 100m 
deep ocean layer about 0.1°K each year (Campbell and Vonder Haar, 
1980). Diagnosing the timing of the hemispheric mean OLWR extremes 
might thus provide significant clues in the long range weather fore- 
cast. 

NOA is estimated by using an empirical formula which relates the 
observed radiance of the 10 y m window channel to the OLWF. A theore- 
tical radiative transfer calculation, rather than observation, is 
used to derive this empirical formula. A recent investigation by 
Ellingson and Ferraro (1983) found that this operational technique 
systematically overestimates OLWF when the technique is applied to 
the radiance from areas with partial cover of high or middle spec- 
trally black cloud. Thus global NOA is about 8 W/m2 higher than 
that of CV. For our study, this systematic over-estimation can possi- 
bly lead to an under estimation of the effect of cloud. A simple 
method to alleviate the difference between CD and the nominal OLWF 



Figure 6.2 Time history of global average OLWR measured by Nimbus 
and 7. (Courtesy of T. D. Bess) 
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is by adjusting GLCLC. However, for the analysis, the temporal and 
spacial variations, as well as CR are not influenced. 

6.2 Zonally Averaged OLWR 

Figures 6.3 and 6.4 show the calculated zonally averaged CO for 
January and July respectively. CR and NOA are also plotted on the 
figures for comparison. 

The typical zonally averaged OLWR pattern, a minimum at the 
tropics and two maxima at the subtropics with negative gradient from 
the maximums toward the poles, are well reproduced. The general 
difference between NOA and CD are within 3.5% depending on the lati- 
tude, except at the antarctic and the minimum at the equator. The 
higher differences in these regions will be futher discussed later. 
The magnitude of the maximum at the subtropics is about 265 to 275 
W/m^, which is about 25 W/m^ higher than that of the minimum at 
the tropics. As discussed in Chapter 4 higher water vapor concen- 
tration at the tropics is one of the major contributor to this mini- 
mum. As well, in the window region, the e-type continuum absorp- 
tion is significant (Warren and Thompson, 1983). From the difference 
between CD and CR, it shows that cloud and RH_ make equiva’,: r '; contri- 
butions to this minimum, which accounts for more than 20W/m2 reduc- 
tion from CR. One of the major cloud regions in this latitude accord- 
ing to Shcrr et_ aj_. (1968) is cloud region 3, which includes 19 out 
of the 36 grids at this latitudinal belt. Table 5.1 reveals that the 
cloud top temperature at this region is much lower than that of the 
adjacent subtropics. The clouds are at higher altitude, 500mb level, 
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compared to 700mb or lower in the neighborhood, while the cloud cover 
changes of the adjacent latitudes are only less than 10%, thus the 
dominant effect in reducing OLWR is the colder cloud top temperature 
and high water vapor concentration. 

At the subtropical maxima, the magnitudes of OLWF are about 
270 W/m 2 with the winter hemisphere slightly higher than that of 
the summer hemisphere, by less than 5 W/m 2 . This is due to stronger 
convection in the summer hemisphere. The difference between NOA and 
CD is smaller than 2.5%, about 5 W/m 2 . Both the OLWR maxima and 
the minima have significant seasonal variations. The locations of 
the maxima, which are associated with the subtropical high pressure 
systems oscillate about 10 to 15 degrees north-south following the 
sun's position. The maximum in the Southern Hemisphere has smaller 
displacement than that of the Northern Hemisphere. The larger ocean 
heat capacity in the Southern Hemisphere stablizes this variation 
(see Figure 6.3). The minimum, which is associated with the ITCZ, 
has less seasonal variation than the maximum. From the observation, 
the ITCZ is essentially located to the north of the equator except 
in January and July. Still, in the CD case, the displacement of 
this minimum is less than that of NOA. This discrepency is mainly 
due to the fixed cloud region classification of Sherr et a_h (1968) 
(see Figure 5.1), which is not seasonal dependent. 

There are larger differences between CD and NOA in the polar 
regions, about 6%, or 6 to 10 W/m 2 depending on the latitude, 
which is much higher than that of mid-low latitudes. Recall from the 
clear sky OLWR discussion in Chapter 4 that in the polar regions the 
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temperature is much lower and the OLWR is much less sensitive to 
the relative humidity. From Figures 4.7 and 4.8 we show that CR com- 
pared best with the other experiments at the polar region, thus the 
difference in the cloudy sky case between CD and NOA is predominantly 
from the effect of cloudiness and cloud top temperature. Since these 
are the places where there are strong surface temperature inversions, 
extending up to higher than 700mb, the effect of inserting a cloud 
deck in the inversion is to increase OLWR. As shown in Figures 6,1 
and 6.2, the difference between CD and NOA increases with latitude, 
which is opposite to the clear sky case. The reason for this feature 
is that to minimize the error in compiling cloud top temperature 
{see section 5.5), each grid value has been weighted with area, such 
that the globally averaged OLWR can be more accurate. In cloud 
Region 24, as well as region 15 at the north pole which extend from 
65°S to 85°S, the weighting of 65°S is cos(65°), about .42, which 
is higher than the sum of the weighting of 75° and 85°, .25 and 
.09 respectively. Thus the cloud top temperature is predominantly 
determined at the latitude of 65°S, so that we have very small 
difference at 65°S between CD and NOA and the difference increases 
toward the poles where there is less weighting. 

Zonal Temporal Variation 

To illustrate time-space variation, a series of time-latitude con- 
tours of T s , RH, CR, cloud cover, CD and NOA are plotted. 

Figure 6.5 shows the time history of zonally averaged climatolo- 
gical T s . The most significant feature is the asymmetrical seasonal 
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variation between the two hemispheres. The strong annual cycle at the 
Northern Hemisphere is due to the large continent distribution which 
contributes rapid heating and cooling, while in the Southern hemi- 
sphere large ocean heat capacity evens out most of the annual cycle. 
In the tropics, the temperature variation is basically non existent 
due to constant solar duration. 

Figure 6.6 shows the time history of zonally averaged ^H. Essen- 
tially RH varies in a very small range, from 0.4 to 0.6 for most of 
the time and space. Higher values in the tropics and lower values 
in the subtropics are well defined; it is evident that the contours 
follow the sun's position. 

Figure 6.7 shows the time history of zonally averaged CR. It is 
clear that CR has a pattern pretty similar to that of the T s except 
in the tropics and subtropics, where the surface temperature is 
sufficiently high. Thus the variation of RW modulates CR in the low 
latitudes. The difference between the maximum in the tropics and 
the minimum in the subtropics is about 10 W/m^. 

Figure 6.8 shows the time history of zonally averaged cloud cover. 
Although there is similarity between the cloud cover and RjH , especi- 
ally in the southern hemisphere, the cloud cover pattern is less 
smooth and spreads over a larger range, 0.4 to 0.8, over the globe. 
The very low value of 0.4 in the Antarctic winter is quite questiona- 
ble. 

Table 6.2 tabulates the time history of zonally averaged T s -T c . 
The corresponding figure is on Appendix A, 17. The tropical regions 
have the highest T s -T c difference, due to strong convective activity. 
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It is also clear that the convective activity follows the position of 
the sun. Although the rate of change of T$-T c is small in tne mid 
latitudes, where most areas are between 10 and 20 °K, the annual 
variation is strong in the Northern Hemisphere and practically non- 
existant in the Southern Hemisphere. Surface temperature inversions 
are constant in Antarctic and polar night time. 

Figure 6.9 shows the time history of zonal ly averaged CD, and NOA 
is plotted on Figure 6.10 for comparison. The model results CD 
agree well with satellite observation NOA, both in magnitudes and 
the contour pattern. Again, basically the feature of CD can be 
decomposed into two parts: in the high latitudes it is the same as 
that of T s and CR, and in the low latitudes it is mciulated by jlH 
and cloud cover. Comparing CD and CR reveals one interesting feature 
that cloud reduces OLWR very homogeneously, about 10 W/m^, for all 
the latitudes but the tropics, where the influence of cloud is much 
stronger, up to 25 W/m?. The asymmetric seasonal variation of the 
two hemisphere is not influenced by introducing cloud. Table 6.4 
shows the global averaged OLWR reduced by increasing GLCLC from .45 
to .60 for July case. As for water vapor, cloud has strongest redu- 
cing effect of OLWR in the tropics. Still only the magnitude of 
OLWR is reduced, the pattern of the contour from CR to CD is not 
significantly changed. 

This feature again stress the importance of the water vapor in 
modulating OLWR, especially in the tropics. Also this suggests 
a new LW parameterization of the type F= a + bT s + c RH as relative 
humidity is highly correlated with cloudiness. A parameterization of 
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this type may improve the performance of a North-type climate model 
in the tropical region. 

6.3 Global Distribution of OLWR 

CD for January and July are shown in Figures 6.11 and 6.12. 
Also, for comparison NOA results for January and July are shown in 
Figures 6.13 and 6.14 respectively. The direction of the discussion 
on the OLWR features will start from the north polar region at the 
zero degree longitude and going west for each latitudinal zone. 

There is good agreement between CD and NOA. The major synoptic 
scale OLWR are well simulated except for a few features, r.amely in 
January the low OLWR center located off the western shore of Mexico 
is located too far north. The main reason for this error is from the 
Sherr et_ a_I_. cloud classification. The cloudiness given by Sherr et 
al . for this region is .28, which is substantially lower than expec- 
ted. In general, this region has stratus cloud resulting from the 
cold north Pacific current. T s -T c ranges from 44°K to 69°K. These 
values are larger than in adjacent areas, which are around 15°K. A 
similar problem occurs at the west coast of South America, where the 
off shore low OLWR center belongs to cloud region 3, see Figure 5.1. 
The cloud cover of this region is .80 which is much higher than 
others in its neighborhood of .45 to .60. Examining the RH field, 
figure 4.9, one would find that it is a low RH region, which is not 
in phase with the cloud cover. Apparently, there is an error due 
to cloudiness. One possibility is that the diurnal cycle of cloud 
cover in this region with maximum stratus in the morning is not well 
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Figure 6.13 NOA map f or JanuJry> ln w/nl 2 
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simulated in the Sherr et £]_. model. We will have some words on the 
diurnal cycle at the end of this Chapter. Another problem occurs at 
the central Pacific low OLWR center, which extends too far east for 
both months. Other than these regions, most of the features are well 
simulated both in magnitude and in location. Since it is beyond the 
scope of this study to rework the cloud region classification, we 
will keep it as is, and discuss the the radiative transfer process 
of the majority of the regions which are well simulated. 

Mid and High Latitude of the Northern Hemisphere 

In January the major low OIWF centers at the high latitudes of 
the Northern Hemisphere are at Greenland and Siberia, with low OLWR 
at about 150 W/m^, which are due to large surface snow-ice cover at 
this season. With these centers a wave pattern is formed along these 
latitudes. A strong ridge is located at Davis Strait and Baffin Bay 
region east of northern Canada. This is from the effect of warmer 
sea surface temperature relative to the adjacent continent, see 
Figure 4.12. The temperature diffference between the ocean and the 
continent is about 10°K. The similar feature with much larger 
scale occurs along the coast of both North America and Eurasia con- 
tinents. Europe has a very mild meridional gradient due to the 
North Atlantic current extending beyond Scandinavia to the Barents 
Sea. Although located at 75°N, OLWR at this region is as strong 
as that of 55°N part of Siberia. In the Northern Pacific, the 
Kuroshio current also enhances the ridge at Alaska and Eastern 
Siberia. Comparison of CD with CR, Figure 4.13, shows that cloud's 


influence on these high latitudes is less than 10 W/m^ and very 
homogeneous. Essentially, the contour pattern does not change when 
clouds are introduced; they only reduce the magnitude by approxima- 
tely 10 W/m^. Since the temperature is fairly low at this region, 
OLWR is insensitive to the moisture. Also, T s -T c is not large, 
being less than 5°K in most regions, and the OLWR is mainly from 
the surface. In July, this wave pattern due to continentality and 
ocean current is completely reversed except for the small section 
from Greenland to Scandinavia. Greenland is still a trough due to 
its icy surface and the Scandinavia area is very zonal and the inten- 
sity is increased to about 230 W/m^. Again, in July, higher than 
40°N the influence of cloud is fairly homogeneous and is about 
10 W/rtr, except at eastern Siberia and Manchuria regions where T $ -T c 
is about 30°K, much larger than in adjacent areas where it is less 
than 1Q°K. This area belongs to cloud region 10 which has very 
high seasonal cloud amount change. Figure 4.14 shows in July, the 
cloud cover is about .70, while at January is about .30. 

Subtropics of the Both Hemispheres 

In January, the strong OLWR in the subtropical high pressure re- 
gions, which coincide the descending leg of the Hadley Cell, can be 
found along the latitude circle of 15°N. The strength of this high 
is more than 270 W/m^. This feature is interrupted in the Phillipines 
region, by the tropical ITCZ in the central Pacific. The OLWR here 
is less than 240 W/m^. In the Southern Hemisphere, there is a counter- 
part of this high located along 30°S, with fairly equivalent strength. 
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This feacure is highly modified by the ocean current and topography. 
On the coasts of South America, the stratus is formed by the cold 
ocean current at both sides of South Africa. This contributes to 
the low OLWR. The high at the west side of South Africa is due to 
the Kalahari Desert. This high can not be seen in CR, but in CD, it 
merely decreases by 5 W/m^ from CR, much smaller than the decreases 
from CR to CD in the adjoining areas, where there is more than 30 W/m^ 
decrease from cold cloud top temperature and higher cloudiness. Over 
the desert, T $ -T c is about 10 & K and A c is about .47, while in 
the southern part of the Congo Basin, T s -T c and A c are more than 30°K 
and .60 respecti vely. As the months change to July, these two sub- 
tropical high belts move northward with the sun. The northern compo- 
nent now covers a much larger area from 20° to about 40°N, depending 
on the longitude. Over the oceans, it is located about 10° further 
south than over the continent. At North Africa, the Sahara Desert 
at low latitudes also enhance its strength to 310 W/m^, which is the 
highest of the world. The southern component is located about 15°S. 
However, its continuity is interrupted by some of the expansion of 
the tropical low OLWR center, such as the Andes of South America, 
eastern Indonesia and the southern part of the Congo. The strength 
is about 270 W/m2, comparable to its northern counterpart. 

Tropics 

At the tropics, the pattern is the least zonal. First, in January 
there are three major low OLWR regions: the Amazon, Indonesia and the 
Congo. Second, Madagascar is a weak low OLWR region. Westward from 
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Central Africa the Amazon has very low OL WR „ about 210 W/m 2 . At the 
southeast Pacific Ocean, it increases to more than 270 W/m 2 then at 
Indonesia, it decreases to about 240w/m 2 . Westward at the Indian 
Ocean monsoonal circulation OLWR increases to more than 275 W/m 2 in 
the clear season. At the African Congo, OLWR decreases to 240 W/m 2 , 
and at the Southern Atlantic Ocean, it increases to more than 260 W/m 2 
again. Here, basically the strength is equivalent to the subtropical 
high OLWR region, which is also true on the southeast Pacific. In 
the Indian Ocean, south of the monsoonal circulation, the OLWR is 
only up to 260 W/m 2 , slightly less than that of the other two large 
oceans. As we have discussed in the sections on the humidity field 
and cloud cover field, the longitudinal variation is strongest in 
the tropical zone. This argument can also be applied to the July 
case, where this region is located north to that described above, 
with some feature changes. It is very clear that in July this low 
belt moves northward following the sun. The India monsoonal region 
is relatively cloudless, and Madagascar becomes a high OLWR center. 
The three lows .'long the tropics move northward. It is not uncommon 
that T s -T c is higher than 55°K in the tropics. For cloud region 3, 
the tropical convective cloud amount for both January and July are 
higher than 60% and has little seasonal variation. The influence 
of ocean current along the continent in the Southern Hemisphere is 
still detectable. Southeast China is under southeast flow f-'orn the 
Pacific high, which brings in ocean moisture. This moist air is 
lifted by the low plateau to create the rain season, and thus a low 
OLWR center. As in the zonal ly averaged case, the minimum of OLWR 
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in ITCZ region stays north of the equator <pt January and February 
when the cloudiness of the three major 1u\ - 5 hi ft to the South. Over 
the open ocean the cloudiness associated with the ITCZ generally 
remains north of the equator. 

At the Aitezon, OlWR decreases from about 267w/m 2 for CR to less 
than 210 W/m 2 for CO. At this region, there is very strong diurnal 
variation with extreme convecti vity. T s -T c at this region can reach 
more than 75°K, with T c less than -35°C. At the peak of this diurnal 
cycle, the T $ -T c can be even larger, T c colder and OLWR below 180 W/m 2 . 
A similar argument can be applied to the Congo, although it is less 
extreme, 265 W/m 2 for CR to 231 W/m 2 for CD. Also the case for Indo- 
nesia is less extreme 273 W/m 2 for CR to 239 W/m 2 . Another feature 
is that the Amazon has more seasonal variation than the other two low 
OLWR centers. It is believed that over Indonesia, the cloud formation 
due to ITCZ has less diurnal cycles, which is not very clear in Sherr's 
cloud model. The ascending leg of the east-west Walker circulation 
coinciding with this convective region causes each cloud cluster to 
last for a few days. This can explain the reason why there is less 
seasonal variation at this region. One very interesting feature is 
that in July, the southern hemispheric winter, there is less convec- 
tivity over the Amazon and the Congo, thus OLWR is higher than that 
of January by about 10 W/m 2 ; over the Congo the difference is even 
larger. Over the strong OLWR region at the tropics, seasonal varia- 
tion is not as prominent as the low OLWR regions; both southeast 
Pacific and central Pacific highs stay fairly constant. 
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Mid and High Latitudes of the Southern Hemisphere 

South of the southern component of subtropical high OLWR band, 
from 40° to 70°S, the feature is very zonal due the fact that only 
ocean is present here. The average gradient of OLWR with the latitude 
is about 1.5 and 1.7 W/m^ per degree latitude for January and July 
respectively, which is more moderate than that of the Northern Hemis- 
phere with 1.8 and 1.1 W/m^ per degree latitude for January and 
July respecti vely. At Antarctic, the OLWR contour is deformed by 
the orography. 

Longitudinal Deviation 

To demonstrate the longitudinal variability, we reconstruct an 
annual longitudinal OLWR deviation contour by calculating: 

Idv (j » k ) ~La n ( J )”£az ) 6.1 

where Fjj v is the annual zonal OLWR deviation for the grid point of 
latitude k, longitude j, £30 is the annual average OLWR, and F_^ z 
is the annual zonally averaged OLWR for the latitude k. This is 
shown on Figure 6.15. One would easily identify that the highest 
deviation are between 5°N and 20°S, with three dipoles of about £ 25 
W/m^. Other regions that have strong deviation are Tibet, Mideast, 
and Greenland. 

6.4 Error Statistics 

The percentage difference between NQA and CD is discretized by 
rounding off to the nearest integer for every grid box. Figures 
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6.16 and 6.17 are scatter plot histograms of the discretized percen- 
tage difference versus numbers of grid boxes for January and July 
respectively. For both cases, the distributions are quite normal 
with standard deviations of 4.6 for January and 4. IX for July. Note 
that the total number of grid boxes is 648. These histograms mani- 
fest the reasonably good performance of this model. Table 6.4 
tabulates the global standard deviation for each cloud type region 
for each month of the model year. Although there is some difference 
between the highest standard deviation, i.e. 4.6X for January and 
February, and the lowest, 3.4X for Septerriber, October and November, 
there is no apparent seasonal dependence. 

To study in detail how the difference is distributed in longitude 
and latitude. Tables 6.5 and 6.6 tabulate January and July percentage 
difference on grid box coordinates accompanied by the Sherr et al_. 
(1968) cloud region classification for easy identification. One can 
find with ease that for most cloud regions there are both positive 
and negative differences in each cloud block. Although these differ- 
ences can be from the deviation of any radiative parameter, the 
larges; error is possibly from the variability of cloudiness. As 
discussed in Chapter 2, our grid box has a size of 10°xl0°. Prac- 
tically, homogeneous cloud can not extend more than a few grid boxes 
of such a size. With mixed clouds in a cloud block, using a single 
cloud top level to specify cloud top temperature is bound to introduce 
error. What is striking is that for most regions, the er^or is 
within 5%. As one would guess, the highest error occurs at the 

tropics with its strong convecti vity, where the errors are higher 

# 


* 



6.4 RMS error in percentage for each cloud type and global standard deviation for each month 
in percentage. Number of grid box for each cloud type is also shown. 
















rr»<-|«(T *f.r fMT* Hit 


144 





NMMK 


Table 6.5 January percentage OLWR difference, (N0A~CD)/N0A, tabulated on 
grid coordinate. JxlO 0 is west longitude, 85°-(k-l) xlO° is 
latitude. Regions are separated according to Sherr et al. 
(1968) cloud region classification tabulated below. 
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Table 6.6 Same a$ table 6.5 but for July 
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than 10%. The strong convectivity in the region generally generates 
non-stratus type cloud, which usually is less homogeneous. One 
example is the cloud region 3 located in the central Pacific and the 
northern part of South America. For the former region, the difference 
ranges from +13% to -14%; for the latter region, the difference 
ranges from +20% to -5%. No other regions have such a diverse dif- 
ference within a cloud block. It is recommended that this region be 
further studied and possibly subdivided and reclassified. One other 
interesting feature is that at the polar regions, cloud region 15 
for the northern one and 24 for the southern one, each one has 
higher error at polar night. This coincides with the error of cloud 
fraction shown in Figures 5.2 and 5.3. Obviously, the cloud observa- 
tion in the polar night time needs to be improved. 

To demonstrate the error according to each cloud type, monthly 
latitudinal weighted RMS errors for each cloud region are calcu- 
lated and tabulated on taHe 6.4. This table 6.4 provides a good 
source for tracing the model performance based on each cloud type, 
and as well indicates the area for future improvement. It is intrigu- 
ing to see that more than half of the cloud types have consistant 
and small, below 5%, errors, while some do not, such as cloud regions 
3, 12, 16, 17, 19 and 25. Except for two of the four grid boxs of 
cloud type 19 that are located at 35°N, all the rest are located in 
the tropics and subtropics region between 25°N and 25°S. This 
again stresses the difficulty of modeling the cloud in the tropics 
region. Another note is that when compared table 6.4 with the cloud 
top level, table 5.7, there exists no correlation. It is recalled 
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from Chapter 5 that cloud regions 3 and 12 have very steady cloud 
top level throughout the year, while cloud regions 17, 19 and 25 
have strong cloud top level variation. 

Diurnal Cycle 

Another area which should be discussed is the bias of satellite 
measurement due to OLWR diurnal cycle. Polar orbiting satellites 
are often used to measure the earth radiation budget. Since the 
satellite passing time is fixed at a certain local time, the diurnal 
cycle must be understood in order to estimate the total budget. 
Unfortunately, the OLWR diurnal cycle changes from region to region. 
Over deserts where there is little cloud, the OLWR essentially follows 
the surface temperature and peaks at mid-afternoon. Over rain forests 
where strong convectivities occur in the afternoon, the cold cloud top 
temperature of cumulus minimizes OLWR. Over the west coast of America, 
where stratus cloud forms in the early morning and dissipates later, 
OLWR has its minimum in the morning. Other types of topography have 
different characteristics. Figure 6.18 from Saunders and Hunt (1980) 
reveals some of the different types of OLWR diurnal cycles. Minnis 
and Harrison (1984), Kandel (1983) and Brook and Minnis (1984) have 
pioneered some of the work in this areas. Still many areas have not 
been investigated. It is important to explore this phenomenon further 
in order to justify the estimation of total radiation budget from polar 
orbiting satellites. 


Longwave flux (W m "*) 
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Local time (h) 

Figure 6.18 The diurnal variation of OLWR for different types of 
surface after Saunder and Hunt (1980). x, over land; 

• , over sea; ♦, over low Cu cloud; a, over high Cu- 
Nb cloud. 
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CHAPTER 7 


CONCLUSIONS AND RECOMMENDATIONS 

A simple and economical earth outgoing longwave radiation climate 
model has been composed for radiation budget study. The model con- 
sists of the the upward radiative transfer parameterization of Thomp- 
son and Warren (1982), the cloud cover model of Sherr et (1968) 
and a monthly average climatology defined by the data from Crutcher 
and Meserve (1971) and Tal jaard et (1969). Additional required 
information is provided by the empirical lOOmb water vapor mixing 
ratio equation of Harries (1976) and mixing ratio interpolation scheme 
of Briegleb and Ramanathan (1982). Reference data from the NOAA scan- 
ning radiometer are used to compile a climatology for the cloud top 
temperature. Both clear sky and cloudy sky cases are calculated 
and discussed for global average, zonal average and global distributed 
cases. 

The clear sky case shows that the OLWR field is highly modulated 
by water vapor, especially in the tropics. The strongest longitudinal 
variation occurs in the tropics. This variation can be mostly ex- 
plained by the strong water vapor gradient. Although in the zonal 
average case the tropics has a minimum in OLWR, the minimum is essen- 
tially contributed by a few very low flux regions, such as the Amazon, 
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Indonesia and the Congo. There are regions in the tropics such that 
their OLWR is as intense as that of the subtropics. In the high 
latitudes, where cold air contains less water vapor, OLWR is basically 
modulated by the surface temperature. Thus, the topographical heat 
capacity becomes a dominant factor in determing the distribution. 

Clouds enhance water vapor modulation of OLWR. Tropical clouds 
have the coldest cloud top temperature. This again increases the 
longitudinal variation in the region. However, in the polar region, 
where a temperature inversion is prominent, the cloud top temperature 
is warmer than the surface. Hence, cloud has the effect of increasing 
OLWR. The implication of this cloud mechanism is that the latitudinal 
gradient of net radiation is thus further increased, and the forcing 
of the general atmospheric circulation is substantially different 
due to the increased additional available energy. 

In che clear sky case, regional results of North and South America 
are compared with GOES data analyzed by Minnis and Harrison (1984) 
with reasonable agreement. The cloudy sky cases also agree well with 
the NOAA satellite observation. T hese cases minutely overestimate the 
global average and have less error in the mid-high latitudes, and 
higher errors in the tropics and the winter polar region in both of 
the zonally averaged and global distributed cases, with standard 
deviation equal to 4.6% and 4.1% for January and July respectively. 

In simulating the case with clouds, the global cloud adjusting 
factor GLCLC defined by Bartman (1980) can be used for minor tuning. 
However, the OLWR fieid can be more effectively adjusted by changing 
cloud top temperature. The current compiled cloud top temperature 
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climatology is probably overestimated due to the inversion scheme 
used for the NOAA scanning radiometer data. Using other satellite 
data as the reference would change the result. 

The climatological cloud region classification by Sherr et aT_. 
(1968) is relevent for simulating the cloud distrbution In most areas 
except the tropical regions where the cloud region classification 
should be further investigated whan new cloud data become available. 
Since many semi-annual air mass systems not only change the intensity, 
but also the area of coverage, seasonality of cloud region possibly 
needs to be included. Why the cloud cover climatology of London 
(1957) and von Loon (1972) significantly differ from that of the 
Sherr et £l_. (1968) in the polar night region should be addressed 
as well. 

This study also presents a global RH climatology. Since the 
cloud cover is highly correlated with the water vapor distribution, 
and the water vapor can be more objectively measured, regardless of 
the fact that there still lack good measurement between 400 and 
lOOmb, some effort should be spent on establishing the empirical 
relationship between the two. Thus, the more simple Budyko-Seller 
type longwave parameterization can be upgraded to the form F=a+bT+cRH 
so as to improve OLWR computation in the tropical region. 

This model can be implemented to study the time-space variation of 
the OLWR field, as well as the sensitivity of the OLWR field to the 
different input parameters. One example of such a case is to overlay 
the current data base with the statistics of the atmospheric variabi- 
lity (Oort, 1983, the data tape would be available soon) as the input 
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data. Variability of the OLWR can thus be analyzed. Another 
application is to study the sensitivity of the earth net radiation 
budget to cloudiness. Since both the Time- Variable Albedo Model of 
Bartman (1980) and this longwave model use the cloud model of Sherr, 
et al_. (1968), this can be achieved first by using the albedo 
model to compute a solar absorbed radiation. The variation of the 
cloud cover to achieve global annual energy balance can then be 
applied to the composite of the two models. 3ecause the absorptivity 
of cloud decreases as the solar zenith angle increases, one would 
again expect that the latitudinal net energy gradient increases due 


to the effect of cloud. 
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R map for Apri I , 
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p for October, in W/m^ 
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A. 10 Cloud cover map for October 
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or April, in W/m 2 
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A. 16 NOA map for October, in v/m 
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A. 18 Comparison of the coefficent values from different 
OLWR parameterizations after Gupta et £l_. (1978); 
see Equation 2.2 for Cj and C 2 , Z c T<F assumed 
black cloud height. 


Reference Cal 

Ci 

cm"2 min” 1 

C 2 

Cal cm”2 min” 1 

Budyko (1969) 
(empirical ) 

0.367 

0.103 

Cess (1974) 
(Z c = 6.8 km) 

0.368 

0.102 

Schneider (1972) 
(Z c = 5.5 km) 

0.399 

0.107 

Gupta et al. (1978) 
(Z c = 6 km) 

0.367 

0.099 

* 100 W/m2 = 0.144 

Cal cm”2 min” 1 
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